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ABSTRACT 

We present Cygnus X in a new multi-wavelength perspective based on an unbiased BLAST survey 
at 250, 350, and 500 /im, combined with rich datasets for this well- studied region. Our primary goal 
is to investigate the early stages of high mass star formation. We have detected 184 compact sources 
in various stages of evolution across all three BLAST bands. From their well-constrained spectral 
energy distributions, we obtain the physical properties mass, surface density, bolometric luminosity, 
and dust temperature. Some of the bright sources reaching 40 K contain well-known compact H II 
regions. We relate these to other sources at earlier stages of evolution via the energetics as deduced 
from their position in the luminosity- mass (L-M) diagram. The BLAST spectral coverage, near the 
peak of the spectral energy distribution of the dust, reveals fainter sources too cool (~ 10 K) to be 
seen by earlier shorter- wavelength surveys like IRAS. We detect thermal emission from infrared dark 
clouds and investigate the phenomenon of cold "starless cores" more generally. Spitzer images of 
these cold sources often show stellar nurseries, but these potential sites for massive star formation 
are "starless" in the sense that to date there is no massive protostar in a vigorous accretion phase. 
We discuss evolution in the context of the L-M diagram. Theory raises some interesting possibilities: 
some cold massive compact sources might never form a cluster containing massive stars; and clusters 
with massive stars might not have an identifiable compact cold massive precursor. 
Subject headings: submillimeter — stars: formation — ISM: clouds — balloons 
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We report on an unbiased survey of the Cygnus X 
(Cyg X) high-mass star formation region, conducted in 
2005 by the Balloon-borne Large Aperture Submillimeter 
Telescope (BLAST), a 2-m stratospheric telesc ope that 



maps sim ultaneously at 250, 350, and 500 /am (Pascale 
|et al.|2008| ). A primary ambition for BLAST was to study 
the earliest stages of massive protostellar evolution. Mas- 
sive stars play an important role in Galactic ecology, ini- 
tially through protostellar outflows and later through ra- 
diation pressure, ionization, stellar winds, and supernova 
explosions, and yet, in spite of their pivotal role, very lit- 



tle is known about their formation ( |Zinnecker fc Yorke 
^007). 



Massive molecular clouds are favorable sites for mas- 
sive star formation. If they are sufficiently dense to 
be self-gravitating, then they also have high extinction. 
Even in the near to mid infrared they have significant op- 
tical depth, affecting observations of any embedded pro- 
tostars, and leading to the " infrared dark cloud" ( IRDC) 
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henomenon seen by ISO dPerault et al. 1996), MSX 



Egan et al.|[l998; Simon et aT|2UU6|), and Spitzer IRAC 



and MIPS fl(Jarey et al.||2UU5[ ~ p^aemer et al.]|2Q10a| |. 
Some IRDCs are nurseries of young protostellar objects 
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and protoclusters (|Simon et al.||2006 ). However, IRDCs, 
by definition found in silhouette, require a luminous dif- 
fuse background. Massive clouds can be discovered di- 
rectly, without this detection bias, by their far-infrared to 
mm-wave optically-thin thermal dust emission. BLAST 
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exploits this fact. 

We introduce the BLAST imaging of Cyg X in 
BLAST detects thermal emission from dust in largely 
neutral regions and so provides a complementary view of 
the interstellar medium (ISM) com pare d to, for example, 
radio emission from ionized gas (§ |2.1| though both are 
ultimately dependent on the local stellar radiation field 
and so are related spatially in some predictable ways). 
An important feature of BLAST imaging is its capabil- 
ity of observing and characterizing extended structures 
in t he s tar- form ing env iron men t at different spatial scales 
(§§|2!3|and|2!4|). In [2^1 and [2l$] we describe the iden- 



tification and quantification of BLAST compact sources. 
Appendix[A]discusses measurement of their flux densities 
at different wavelengths for use in the multi-wavelength 
spectral energy distribution (SED). 

Blind submillimeter surveys w ith BLAST ( Chapin 
et al.||2008| |Netterfield et al.||2QQ9|>, precursors to those 



with Herschel (e.g., Molinari et al. 2Q1Q[), are ideally 



suited for finding and characterizing cold sources, be 
cause of the designed coverage near the peak of the cold 
dust spectral energy distribution (SED). This allows us 
to determine where the dust temperature T is low, a 
key requirement before determining the cold dust col- 
umn densities and masses. BLAST of course also sees the 
later evolutionary stages when the dust is warmed up by 
the forming protostars, and again the mult i- wavelength 
coverage tightly constrains both dust temperature and 
bolometric luminosity (Lboi)- 

In order to interpret the BLAST emission more fully, 
and place the compact sources in context, we make use of 
many different surveys of Cyg X in other tracers, includ- 
ing both continuum emission and molecular lines. These 
rich mult i- wavelength, multi-species surveys of Cyg X are 
the basis of a brief overview of the physical environment 
in§! 

The spectral energy distributions (SEDs) and deduced 
properties T, L, and mass M of the com pact sources are 
presented in § |4j As we discuss in § |4.6[ the mass of the 
compact source, along with its luminosity (which maps 
into T), determine its position on an evolu tionary track 
in the L-M diagram ( [Molinari et al.|[2QQ8|) . Investigat- 
ing the stages of high mass protostellar and protocluster 
evolution is a prime aspect of this paper. 

In the earliest stages, before there is any (significant) 
internal source of energy from accretion or nuclear burn- 
ing, the dust is heated only by the external radiation 
field. Because of the high extinction the dust is cold, 
with T typically ~ 15 K or lower, but in any case colder 
than the "ambient" temperature of diffuse dust in the 
same external radiative environment (which near these 
star forming molecular clouds might be much more in- 
tense than the local interstellar radiation field). Such 
cold clouds, and the initial gravitationally-bound con- 
densations within them are therefore rather invisible to 
IRAS, but are in principle detectable by BLAST and 
grou nd-based surveys targeting regions of high extinction 
(e.g., |Motte et al. [20071 using MAMBO at 1.2 mm). Ac- 
cording to the time sequence and nomenclature used for 
low mass star formation, sources detected in this earlier 
evolutionary stage could be called "class —1." However, 
as a reminder of the energetics, we prefer the physically- 
motivated shorthand "stage E" ( from Externally heated, 
but also usefully Earliest - see § |5.6[). Likewise, we call 



sources in the subsequent evolutionary stage, when there 
is sufficiently vigorous accretion power internally to raise 
the d ust to the ambient or higher temperature, "stage A" 
(see § 5.5 the low mass analog would be "class 0"). 

An illustrated view of the evolutionary stages of mas- 
sive star formation, illuminated from the BLAST sub- 
millimeter perspective, is presented in § |5j 

A key factor that is important to the interpretation 
of our particular observations is th at massive stars form 
in stellar groups and star clusters ( |Lada fc Ta da 2003). 
The approximately V angular resolution (EWHIVlj ol the 
BLAST images corresponds to 0.5 pc at 1.7 kp c. This 
is the size scale (sometimes called a "clump" - Bergin 



fc Ta falla 2007]) of the obser ved embedded and newly- 
emerg in g clu sters in Cyg X ( |Le Duigou fc Knodlseder] 
|2002[ § |3.1[) . Thus it is relatively straightforward to 
search for trie precursors to these protoclusters, within 
which massive protostars will be the dominant source 
of luminosity and ionizing radiation. However, resolving 
protoclusters into unconfused condensations (a smaller 
linear scale sometimes called "cores") that are destined 
to form single (proto) stars requires more angular resolu- 
tion than has been available in surveys (even with the 
11" beam of MAMB02 which is smaller than achieved 
here with BLAST). Submillimeter observations with the 
larger telescope of Herschel will be better (18" at 250 jam 
with the SPIRE camera), but not immune to these con- 
siderations. 

Empirically, |Bontemps et al.| ( |2009| ) have made an ex- 
tensive study ol fragmentation and sub- fragmentation in- 
side targeted molecular dense cores in the Cyg X region 
with interferometric observations at Plateau de Bure. 
Theoretically, it is still unclear how clusters containing 
massive stars form ([McKee fc Tan|2003l|Krumholz fc Mc- 
| Kee||2008| [Bonnell et al.||i^7[ |^L)0ij |20Q7[ |&mith et al. 
2009 ; [Wang et al.' 2010p . Some insight provided by the 
BLAST observations and ancillary data, coupled with 
current theory is the focus of the discussion in § [6j 

2. BLAST IMAGING OF Cyg X 

Cyg X, positioned in the Galactic Plane at about 
/ = 80° an d about 1.7 kpc away from the Sun ( |Schneider 
|et al.|2006| ) , has long be en known for its massive star for- 
mation. With BLAST ( |Pascale et al. 
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2008) we surveyed 10 deg 2 in Cyg X for 10.6 hr dur- 



ing the June 2005 flight (BLAST05), mapping the area 
on three visits to provide cross-linked scanning. In ad- 
dition, a significant amount of calibration-related time 
was spent observing a circular cap of radius 1° centred 
on W75N, resulting in even higher signal to noise and 
dynamic range there. 

Final maps on 15" p ixels were produced u sing the 
map-maker SANEPIC flPatanchon et al.| [2008] ). The 
combination of high scan speed and low 1/ f knee, to- 
gether with the cross-linking and common-mode removal 
in SANEPIC, produces a map retaining diffuse low spa- 
tial frequency emission, in contrast to current ground- 
based mm- wave mapping which contends with the atmo- 
sphere and therefore is spatially filtered and emphasizes 
the compact structures. However, preprocessing of the 
time-ordered data to remove very low frequency drifts 
makes the SANEPIC map average near zero, the DC 
level having been effectively removed. 

BLAST05 was designed to produce diffraction-limited 
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Fig. 1. — Upper: enlarged section of 250 jum BLAST05 image 
in Fig. [2] Lower: deconvolved image using a Lucy- Richardson 
algorithm. Full width half maxima of the deconvolved sources are 
about V . 



beams with FWHM 30", 42", and 60" at 250, 350, and 
500 /im, respectively. However, BLAST05 suffered from 
an uncharacterized failure relating to the primary mirror 
figure and focus. Nevertheless, even with the corrupted 
beam the resolution is about 3' full width half power 
flTruch et al.|2008||Roy et al.|2010[ ), comparable to IRAS, 
though disappointing compared to th e diffraction-limited 



images obtained in th e 2006 flight ( |Truch et al.||2QQ9 
Netterfield et al.|2QQ9| >. 



The characteristic beam shape can be seen clearly at 
many locations throughout the field, indicating the pres- 
ence of point sources with angular size comparable to the 
diffraction limit. Obviously we want to work with images 
with as good resolution as possible. We removed the 
worst effects of the corrupted beam (see Roy et al. 2010, 
in prepar ation) by a Lucy- Richardson deconvolut ion al- 
gorithm (Richardson 1972 Lucy|[l974| ). The significant 
improvement is illustrated in Figure [T] Even so the effec- 
tive resolution is about 1', less dependent on wavelength 
than for diffraction (which has the silver-lining benefit 
that beam matching is not an issue). Moreover, there 
are artefacts introduced near very strong sources. Math- 
ematically, Lucy-Richardson deconvolution should con- 
serve flux, and we have verified this through simula tion s 
and measurements on isolated compact sources (§ 2.6). 
The astrometric a ccura cy is good, about 5" when the de- 
convolved maps (§ 2.6) are aligned with compact sources 
in comm on with M7TMBQ2 jMotte et al.||2QQ7| > and/or 
SCUBA pi Francesco et al.l2DD5| r" 

The full 250-/mi map is shown in the upper panel of 
Figure [2] BLAST measures the dust continuum emission 
along the line of sight and probes the major star-forming 
sites in the region by detecting reprocessed stellar or pro- 
tostellar radiation. The positions of luminous s tars in 
the Cyg OB association (iNegueruela et al.' 2008) and a 



ence. Expanding H II regions, bubbles, and stellar out- 
flows, and the interplay with structures in the interstel- 
lar medium make the morphology of extended emission 
in the Cyg X region quite complex. A comprehensive 
review coveri ng important physical aspect s of Cyg X is 
presented by |Reipurth fc Schneider| ( |2008[ ) . 

2.1. A Complementary View using Radio Continuum 
Emission 

A complementary way to see the effects of the high- 
energy radiation from hot massive stars is via radio emis- 
sion from the ionized gas. The lower panel of Fig ure [2 
shows the 2 1-cm continuum map from the CGPS (Taylor 
et al.|200"3| ), which has V resolution, about the same res- 
olution as the BLAST maps. A number of H II regions 
are l abelled for reference : num erals n are DR n regions 
from Downes fc Rinehart] ( |1966[ ), ECX6-* from I Wendker| 
al. ( |l99l| ), and otherwise G/ (accurate coordinates for 
these three] regions being G79.957+0.866, G80.451+1.041 
and G81.445 + 0.480). Note that the recombination line 
velocity for G79.9 is -14.8 k m s" 1 and for ECX6-27 is 



-64.4 km s -1 |Lockman||l989 ) , putting these beyond the 
Cyg OB2 association. 

In diffuse regions the free- free emission is optically thin, 
but for the most dense compact H II regions, it can be 
self-absorbed. There is rough correspondence with the 
dust emission, as might be expected given the impact of 
OB star radiation on both. Note, however, that unlike 
the dust emission, which depends on the column den- 
sity, the H II emission depends on the emission mea- 
sure, involving the square of the local density. Further- 
more, ionizing radiation is obviously essential to produce 
the plasma. This radiation, the rest of the non-ionizing 
starlight, trapped Lyman-a, and cooling lines all heat 
any dust surviving in the plasma. Most of the sub- 
ionizing radiation escapes the plasma to heat any dust 
outside the H II region. In an edge-on geometry, as oc- 
curs in extended regions like DR22, the warm dust imme- 
diately outside the arc-shaped ionization front (I-front) 
can be seen (Fig. [2] and magnified views in figures be- 
low). In a related way, in DR23 the ionization is clearly 
stopped on the right by a dense cloud which forms a 
"bay" in the radio emission, and a corresponding bright 
feature in the dust map. 

There are also large sc ale ionized structures denoted 
CXRn (Cyg X Ridges) by | Wendker eTld] ( p9fl . Those 
in the BLAST region are labelled in Figure [2] These 
are seen dimly in the BLAST image because they have 
relatively low column density (they also do not show very 



number of stellar clusters (§~ XT) are indicated tor refer- 



strongly in the extinc tion map; § 3.4). A much better 
dust emission tracer (§ |2.3| ) of these ionized ridges, indeed 
of all of the diffuse ionized gas, is the 24 /im MIPS image , 
from the Spitzer Cyg X Legacy Survey ( HometaLj2009 ) , 
which looks remarkably like the radio continuum image. 

2.2. Planck Cold Core Survey 

In the midst of this turmoil, we search for relatively 
cold structures which could be the precursors of the next 
generation of stars. This is also th e goal of the lower - 
resolution Planck Cold Core survey flJuvela et aL|2010| ), 
an unbiased search using the Planck all-sky maps of sub- 
millimeter dust emission. Cold cores are characterized 
by a lack of corresponding IRAS 100 /im emission. Be- 
cause the Planck high-frequency angular resolution is 4', 
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Fig. 2. — Upper panel: BLAST05 deconvolved 250 /im dust emission map of Cyg X. Contours are of 21-cm radio c ontinuum emission 
from the CGPS in the lower panel. The large dashed circ le of radius 60' shows th e extent of the Cyg OB2 association ( |Knodlseder|2 000 ) 
and the stars are t he m ost luminous O stars in the list of Negu eruela et al.] ( |2008| > (see text). Circles denote the positions (not extent) of 
star clusters (see § |3. 1| for nomenclature). Lowe r panel: The radio image m turn has contours from the BLAST05 image. Prominent H II 
regions and diffuse ridges are marked (see § |2.1|for nomenclature). 
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clumps - even if more extended than the BLAST clumps 
- will appear as point sources, and smaller structures will 
be beam-diluted. Even for relatively nearby regions like 
Cyg X, the Planck beam would correspond to a spatial 
extent of about 2 pc. 

The highest-frequency Planck channels correspond ap- 
proximately to 350 and 500 jam. Prior to the availability 
of Planck imaging data, we explored the possibilities by 
convolving the two longer-wavelength BLAST images to 
the Planck resolution and combining, them with IRAS 
100 jam in a three-color image (Fig. [3]). At this resolu- 
tion, BLAST stage E sources like C1I6 (§p]and § |5.6.1 ) 
still stand out as being cool. More extended, elo ngated 
structures like the cool IRDC ridge nea r DR1 5 (§|5 .6.2) 
and the DR17 molecular pillar ( 



3.6 



and § 5.1.ip are 



also preserved. Of course the regions of warmer dust are 
highlighted in this image as well. 

2.3. Diffuse Emission 

BLAST maps measure surface brightness (MJy sr _1 ), 
hence column density of the dust M&\ 



I v = M d n v B v (T) = r/am u N R K„B„(T), 



(1) 



where r is the dust to gas ratio, Nn is the H column den- 
sity, and k v is the mass absorption coefficient of the dust. 
For later characterization of the SED we note that the 
spectral dependence of k v is usually scaled with respect 
to a fiducial frequency as k u /ko = (v/vo)^- Note also 
that for absolute measures of column density we would 
have to restore the zero point (DC level) of the map s, as 
we did for the Cas A region ( Sibthorpe et al.||2010[ ) but 
not here. 

In the Galactic Plane, cirrus-like structures in the form 
of bright diffuse emission , termed "interstellar froth" 
( Waller fc Boulanger||1994 ), are distributed everywhere. 
With the com pact structures removed, eith er explicitly 
or by clipping ( Miville-Deschenes et al.|2007 ), the diffuse 
emission has a cirrus- like power spectrum, with lower 
fluctuations at high spatial frequencies. In fact, we have 
found that for this region the power spectrum is quite like 
that for IRAS 100 /im, with amplitude simply scale d by 
the SEP 2 for the appropriate cirrus temperature (Roy 
eTaLl[2QT0 ) 



Somehow, through turbulent motions, thermal insta- 
bilities, and phase transitions, massive localized regions 
become self-gravitating and so at high spatial frequen- 
cies they stand out from the dwindling cirrus structure 
as compact sources. A distribution of masses is expected, 
but the cirrus confusi on l imits our ability to distinguish 
lower mass sources (§ |2.5| ). This situation improves with 
angular res olution, but even fo r Herschel it is a dominant 
limitation flMartin et aL|2010| ). 

Actually, the diffuse emission in the submillimeter and 
mid-infrared wavelengths comes from three different dust 
compone nts, distinguished principally by thei r size distri- 



bution (Desert et al |1990l|Li & Draine|2001l Compiegne 



et al., in preparation): Big Grains (h^Gs), Very Small 
Grains (VSGs), and polycyclic aromatic hydrocarbons 
(PAHs). The BGs, which account for the most of the 
dust mass and therefore most of the longer wavelength 
emission, are in thermal equilibrium. The VSGs have a 
relatively lower share of the total dust mass, even smaller 
in dense regions. They are small enough to experience 



non-equilibrium heating and so broaden the spectrum 
toward shorter wavelengths, beyond the spectral peak of 
the BG emission. Compared to what is expected from 
equilibrium BGs alone, this excess non-equilibrium emis- 
sion appears typically at 60 /im and shorter wavelengths. 
PAHs are also non-equilibrium emitters, in strong spec- 
tral bands pumped by far-ultraviolet radiation. The rel- 
ative amounts of emission i n these three component s can 
vary from region to region. Schneider et al. (2006) used 
the MSX band A images (8.3 /xm; Carey et al. 2005]) ef- 
fectively as a tracer of PAHs and the envi ronment. We 
use S pitzer IRAC band 4 images (8.6 /im; Wer ner et al.| 
2004 ) from the Cyg X Lega cy Survey ( |Hora et al.||2009[ ) 
tor the same purpose (§ 3.3), since they have better sensi- 



tivity and resolution. Spitzer MIPS 24 /im images trace 
VSG emission generally, and the hotter grains associated 
with ionized gas, and so, as commented above, look re- 
markably like the radio continuum images. 

2.4. Dust Temperature from the Diffuse Emission 

From the above equation it is clear that the BG emis- 
sion that BLAST sees is modulated by the dust temper- 
ature T which is in equilibrium with the local radiation 
field, whether in the diffuse medium or within a com- 
pact source. We concentrate here on the former. For our 
exploratory work here on the BGs, we adopt a single- 
temperature SED and therefore fit only data at 100 /im 
and larger to avoid contamination by VSG emission. 

For a sufficiently large and homogeneous region an es- 
timate of the characteristic temperature can be obtained 
by fitting an SED to the square root of the amplitudes of 
power spectra. By this method u sing data at 100, 250, 
350, and 500 /im, Roy et al. (2010) obtained a cirrus dust 
temperature of 19.9±1.3 K for a relatively uncrowded re- 
gion (basically the right half of the BLAST map). How- 
ever, this more global power spectrum method is not 
practical when the region considered is relatively small, 
inhomogeneous, and/or focusing primarily on particular 
sub-structures. 

An alternative approach for smaller regions begins with 
pixel-by-pixel correlations of images with respect to some 
reference image (here 250 /im). Small scale structures 
are remarkably well correlated across the three BLAST 
bands. The slopes of these correlations describe the rel- 
ative SED of the spatially-varying dust emission that is 
changing in common across these images. Note that this 
approach cannot give the dust temperature pixel by pixel 
(that would require the DC offsets too). On the other 
hand, the emission in each pixel is from dust of different 
temperatures along the line of sight, whereas our ap- 
proach isolates and characterizes certain more localized 
spatial components. 

For quantification of the size of the changes in temper- 
ature, we selected regions that appeared to have different 
colors in a map like Figure [3J 

The first region, shown in 500 jam emission in the left 
panel of Figure [4j is within the relative void created by 
the Cyg OB2 stars, near CXR12 (Fig. [2|. Here there 
is relatively strong 100 and 60 /am emission. When the 
SED (left panel of Fig. |5| is fit using j3 = 1.5, which 
is consistent with the value adopted for the analyses of 
compact sources in this paper, the derived dust equilib- 
rium temperature T = 29.7 ± 2.3 K. For /3 = 2, appro- 
priate for local diffuse dust in the atomic gas at high 
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Fig. 3. — Three color BLAST plus IRAS image of Cyg X (500, 350, and 100 /im are represented by red, green, and blue, respectively). 
Images are first convolved to the Planck high frequency 4 ; resolution. 

swept up. However, the geometry is unclear, not im- 
mediately recalling a classical edge-on photodissociation 
region (PDR) in a molecular cloud. 

The selected cooler region, judged from relatively low 
100 and 60 jam emission, is shown by a circle in the right 
panel of Figure [4] Its SED (right panel of Fig. [5]) yields 
an equilibrium temperature of 21.6 ± 0.5 K. Thelower T 
would imply a lower effective interstellar radiation field, 
due to the high extinction, with Ay for this molecular 
structure (§|3 5[) up to 6 mag above the local background. 
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Fig. 4. — Two small regions in Cyg X selected on the basis of a 
three-color 500, 250, and 100 fim image to have warm (left) and 
relatively cold dust (images here at 500 /mi). Note that as in other 
figures, white represents lower intensity. 



latitudes ( jBoulanger et al.||1996| |Li fc Draine||2QQl[ ), T 
is 25.4±2.3 K. This can be compared to the local high- 
latitude equilibrium temperature, 17.5 K. The higher T 
in Cyg X results from the much higher effective interstel- 
lar radiation field from the OB association. The ridge of 
warm dust emission and this particular selected warm 
region above and to the left are also prominent in dif- 
fuse H II emission (Fig. [2]) and 24 /im emission, which 
supports the view that the local radiation field is high 
he re. This region is particularly devoid of CO emission 
(§ |3.5| ). However, the material emitting in the submil- 
limeter must be neutral, because the associated extinc- 
tion AAy is up to 3 mag above the local background 
(and i s spa tially well correlated with the BLAST emis- 



2.5. Compact Sources 

When I v in equation is integrated over a suitable 
solid angle, accounting for the background, the flux den- 
sity S v (typically in Jy) of a "source" is obtained. Like- 
wise, the source mass M (gas plus dust) is related to the 
integrated column density when the distance D is known, 
so that 



S v = L u /(4irD 2 ) = MD~ 2 rn v By{T). 



(2) 



sion; 



3.4p . To the lower right of the CXR12 ridge, the 
and there is some CO emission. 



dust is cooler, and there is some CO emission. CXR12 
is oriented roughly tangentially to the direction to the 
center of the OB association, and might involve material 



Note that "source" refers to the dust emission that is 
being observed, not what heats the dust, and so in par- 
ticular does not imply that there is an embedded source 
of energy (a star or protostar) within the volume of dust 
being measured. 

The definition and characterization of what is a com- 
pact source is a difficult challen ge, and different ex - 
traction sche mes, like Clumpfind I Williams et al. 1 994|), 
Gaussclump flStutzki fc Gues ten 1990; Kramer e t al. | 
1998), or the multi-scale method ol |lVlotte et al.| ( |2UU7[ ), 
can produce different catalogs and source characteris- 
tics. Operationally what is called a "compact source" 
is often basically a structure comparable in size to the 
beam. This immediately cautions that what is a com- 
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Fig. 5. — Left: Relative SED normalized at 250 /im corresponding to the circular region shown in the left-hand image of Fig.E] Triangle 
shows the relative value unity at 250 /im. The SED fitted using (3 = 1.5 has a temperature of 29.7 ± 2.3 K. The 60 /im relativebrightness 
(represented by a circle) is not used in the fit. Right: Relative SED corresponding to circular region shown in right-hand image of Fig. [J] 
The fit temperature is 21.6 ± 0.5 K. 



pact source to one instrument could become resolved into 
multiple components at higher resolution, as has been il- 
lustrated by submm/mm interferometric observations of 
IRAS protoclusters and ho t cores (e.g., with the SMA 
and IRAM Platea u de Bure, |Beuther et al.|2007a|b| [Bom 
temps et aL|2009| ). 

As mentioned above, massive stars form in clusters. 
What BLAST can hope to detect in this survey at V res- 
olution is unlikely to be a single protostar or its progen- 
itor, but rather a protocluster or stellar nursery possibly 
containing many unresolved objects. Of course, depend- 
ing on the mass function, a single object among these 
might dominate the luminosity and ionizing flux. 

We have taken a pragmatic view that since whatever 
underlying structure that is smeared out as a compact 
structure does appear quite like a Gaussian, this is an 
adequate model for extracting flux densities. Further- 
more, the deconvolved beam for this BLAST05 survey is 
of similar size at all three bands, so that the same volume 
of dust is being examined. A corollary is that when we 
use other initially higher-resolution data to expand the 
spectral coverage of multi-wavelength photometry (Ap- 
pendix [A]), we measured the flux density consistently by 
first convolving these maps to the BLAST05 resolution. 

2.6. Compact Source Identification and Photometry 

For compact source identification we used the IDL- 
implemented DAOPHOT "FIND" routine, keeping as 
candidates detected peaks that were above the local 
background in both the 250 and 350 jam images by 3<r 
(empirical local rms fluctuation) . Both source identifica- 
tion and photometry on individual sources are affected 
by the presence of cirrus noise; for example, at 250 /im 



the 1 — a cirrus noise level was evalua ted to be ~ 4 Jy 



in the fainter right portion of the map ( |Roy et al.||2010[ ) 
although it is expected to vary over the map according to 
the median brightness of the region (making catalog com- 
pleteness a function of position in the map). The 500 /im 
map was not productive for further independent identi- 
fications due to at least as great apparent relative confu- 
sion noise. In fact, for 14 faint sources among the above 
it was not possible to measure a reliable 500 jam flux 
density. Given this noise and also the deconvolution arte- 
facts near the brightest sources, we visually checked our 
candidates from the automated list. Since we were not 
depending on a complete survey - for example, to find 
a mass function or relative lifetimes of different stages 
- we were conservative in rejecting sources and a few 
were even added to the list. We were most concerned 
with having sources with accurately determined SEDs, 
so that reliable physical parameters could be determined 

(§B- 

.^otometry was carried out using a Gaussian non- 
linear fit, with amplitude, FWHM, position, and lin- 
ear background as free parameters. Model-independent 
aperture photometry for isolated sources gave the same 
fluxes within 7%, confirming that this is an acceptable 
approach. Cyg X has a sufficiently flat background lo- 
cally that fitting a linear background is an adequate 
approximation. However, in a complex crowded star- 
forming region like much of Cyg X, it is often not possi- 
ble to extract fluxes by aperture photometry or by fitting 
a single Gaussian. In such cases, we adopted a multiple 
Gaussian technique, fitting simultaneously to all candi- 
dates within a 2.5' radius about the parent source. If a 
parent source in a crowded region had a low signal to 
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noise ratio compared to its neighbors then we fixed their 
positions and/or FWHMs as required. 

Positions and flux densities for 184 sources are given 
in Table [T] These Cyg X sources are referred to below 
as 'Cn' where 'n' is an up-to-three digit number, ranging 
from to 183. 

3. CONTEXT: GLOBAL ENVIRONMENT AND 
STRUCTURAL DETAIL 

3.1. Cyg OB 2 Association and Star Clusters 

Cyg OB2 is among the more impressive assemblies of 
OB stars in the Galaxy. It is more compact than a nor- 
mal association and rich enough to have been called a 
"young globular cluster" (]Kn6dlseder||2000|) . The latest 



census indicates that it c ontains some 60-70 O-type stars 
QNegueruela et al. | |2008[). From their semi-observational 
HK diagram, [N egueruels 



et al 



(2008) 



IK diagram, |lNegr 

1.5 kpc, with 1.8 kpc being noticeab. 
adopted dis tanc e of 1.7 kpc, consistent with Schneider 



favor a distance 
y too l arge . Our 



et at] ( 12006) and |Motte et al.| ( |2007) , is close to this and 
a 107o uncertainty in distance is of no consequence in the 
analysis that follows. 

The most luminous hot stars (selected as Mb Q i < —9.2, 
logT eff > 4.56 and M bol < -9.9, 4.46 < log T eff < 4.53) 
are plotted in Figure |2] with position s from |Massey "fc 
Thompson! (|1991D andJComeron et al.| (|2002|). Note that 



star A37 (at G80. 240+0. 135) has what appears to be a 
bow-shock to its right in the Spitzer image s and so might 
be a r unaway star from LK8 in DR15. Kobulni cky et al. 
( |2010| ) studied the bow-shock morphology of A37 (and 
others) , placing this star slightly behind Cyg OB2 at a 
distance of 2.1 kpc 



According to Knodlseder (2000), Cyg OB2 extends to 
a radius of 1° ('60 pc), with half-light radius 13' (6.4 pc). 
Examination of the sub-structure reveals two open clus- 
ters in the core (BBD1 and BBD2; |Bica et al.|2003| ), each 
of radius ~ 2' yet separated by only 6' and appearing to 
form a physical pair. 

A number of highly reddened OB clusters and stel- 
lar groups have been found in the surround ing molecular 
clouds using 2MASS. Dutra fc Blca| ( |200lD performed a 
targeted search toward large angular size H II regions in 
Cyg X; eight clusters and four stellar groups are in this 
BLAST survey region. Five of these have been re c overed 
and characterized by Le Dui gou fc Knodlseder| ([2002), 
including the half-population radius R^o and ce ntral den- 
sity. They als o describe two new clusters. |Comeron| 
Torra ( 199 9]) studied D R18, finding a stellar group. 



and Comeron ^et al. (2002) targeted other ECX (Wend- 
ker et al.|l99l[ ) compact H II regions finding four clusters 



within the BLAST area. The new one, ECX6- 27, has a 
recom bination line velocity of —62.4 krns" 1 (Lockman 
1989) and so is much more distant than Cyg OB 2, more 
like 10 kpc. It has a "double concentration" of stars 



(Comeron et al. 2002[ ) and also a complex d ouble struc- 
ture in the Spitzer images. The catalog of |Bica et al 
([2003] ) also includes the W75N and DR21 (W75S) in- 
frared clusters. In total there are thirteen clusters and 
five stellar groups (see Tablej2|. Their positions (not ex- 
tents) are shown in Figure [zpwith names in precedence 
LK, DB, and ECX6-* where the lists overlap. In sub- 
sequent figures, to indicate the cluster angular size we 
use i? 5Q for the OB star po pulation where available ( |Le| 



3.2. The Radiative Environment: Ionization 

The signature of radio emission in this region of the 
Galactic Plane, as distinct from the radio galaxy Cygnu s 
A, led to the name Cyg X QPiddington fc Minnettfl 952). 
A more contemporary large scale view is given in Fig- 
ure 5 of Uyaniker et al. (2001 ). Radio surveys conducted 
at different frequencies have catalogued the ma in H II 
region complexes (1390 M Hz, | Wester hout||1958| 5 GHz, 
Downes fc Rineha rt|1966| multilrequency, [Wendker et al.| 
1991). What we have observed with BLAST is only the 



central portion, all projected within the Cygnus super- 
bubble. The corresponding 1420 MHz image from the 
CGPS, including the diffuse emission, is shown in Fig- 
ure n 

The free-free emission traces the compact sources and 
extended structures like I- fronts, elephant trunks and 
pillars, and material that is being sculpted by winds. 
Figure [2] shows a striking ridge of emission CXR11 ex- 
tending downward from DR15. Sculpting seen in the 
Spitzer images indicates illumination from the upper left 
(Cyg OB2), but like CXR12, this does not appear to be 
a classical I- fro nt, because the expected neutral tracers 
in the PDR (§ 3.3) are not present immediately to the 



Duigou & Knodlseder 2002), otherwise a radius of V . 



right. 

In a closed geometry, the bolometric luminosity of the 
dust emission is a calorimeter for radiation from any em- 
bedded (proto)stars. Likewise, the spatially-integrated 
radio flux from the same region is a measure of the num- 
ber of ionizing photons. Thus one can check for consis- 
tency, which is particularly useful in assessing the distri- 
bution of masses within an unresolved cluster. For ex- 
ample, dividing the luminosity among several lower mass 
stars will produce a lower ionizing flux. There are com- 
plications of course, in addition to the relative covering 
factors. For example, one must know the distance. If the 
distance is overestimated, then the luminosity will imply 
more massive stars and so a relatively higher ionizing 
flux ( jComeron fc Torra |2001[ ). Furthermore, the radio 
emission is also affected by self-absorption. Thus the 
ionized zone around young OB stars deeply embedded 
in dense gas (ultracompact and hypercompact H II re- 
gions) can appear to be underluminous. Multifrequency 
observations (not pursued here) can of course reveal self- 
absorption by its spectral signature in the radi o SED (see 
the discussion of DR21 in | Wendker et al.|1991[ ), or by the 
high brightness temperature if the source is resolved. 

3.3. The Radiative Environment: Far Ultraviolet 

The non-ionizing radiation heats dust beyond any ion- 
ization front (I-front). This is the main source of the 
extended submillimeter emission surrounding the H II 
regions. As described above, in cases with "edge-on" ge- 
ometries, like in DR22, there is a clear displacement of 
the BLAST emission to the side of the I-front away from 
the ionizing stars. 

The non-ionizing radiation induces non-equilibrium 
emission when the relevant species, VSGs and PAHs, are 
present. Diffuse PAH emission is well traced by MSX 
band A; the MS X map of Cyg X is discussed by |SchneT| 
der et al. (2006). PAH emission can now be seen in more 
exquisite detail in the IRAC band 4 images. Where the 
geometry is favorable, this reveals impressive PAH-fronts 
(defined by where the exciting radiation is attenuated) 
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outside the I-fronts. There is also a sweeping arc to the 
lower right of DR15. This is not immediately adjacent to 
the ionized ridge CXR11 mentioned above (thus not the 
standard I-front-PDR geometry). The arc (and ridge) 



can b e seen as well in MIPS 24 /im emission (iHora et al 



|2009| ) , which we interpret there as VSG emission" 

MIPS 24 jam dust emission also turns out to be a good 
tracer of the plasma, although because of the different 
physics there is not a complete morphological or bright- 
ness correspon denc e. In the figures for the selected re- 
gions below (§ 3.6), we show the 24 /im emission over- 



plotted with 21-cm radio emission contours to illustrate 
their correlation, and also the relationship to clusters and 
the BLAST compact sources. 

3.4. Extinction 

Another tracer of column density is near-infrared ex- 
tinction (often expressed as its equivalent in Ay), which 
can be estimated from the colors of 2MASS sources. Such 
a map by S. Bonte mps (private comm unication) is pre- 
sented in Figure 1 in |Motte et al7] ( |2QQ7| ). A v and BLAST 
emission are well correlated. This correlation can be used 
to calibrate the dust op acity (Martin et al. in prepara- 
tion). In their Figure l, |Schneider et al.j (|2006) show an 
Ha image which strikingly shows the foreground optical 
extinction called the Great Cygnus Rift. 

3.5. The Molecular Reservoir 

Numerous surveys of the Cyg X giant molecular cloud 
(GMC) have been carried out in the molecular line tran- 
sition 12 CO J = 1 -> (Cong 19 771 iDame et al.||l987 
|Leung fc Thaddeus|[l992l ). |sdmeider et al.| ( | 2 (J Un report 
on an extensive multi-transition survey with KOSMA 
( 13 CO J = 2 -> 1, 3 -> 2, and 12 CO J = 3 -> 2) to study 
in detail the spatial structural variations and physical 
conditions. The relative intensity of 13 CO J = 3 — >• 2 to 
13 CO J = 2 — > 1 depends on the local density and ki- 
netic temperature. Their 13 COJ = 2^1 channel maps 
show elongated filamentary structure at scales of 10' to 
20', and smaller sub-structures (CO clumps) embedded 
in the larger cloud fragments. The peaks of the CO emis- 
sion profiles occur over a wide range of radial velocities 
surveyed from +20 to —20 km s _1 (there is not much gas 
beyond —10 km s _1 and gas at more negative "Perseus 
a rm velocities" is not in cluded) . 



Schneider et al.| ([2006]) provide a comprehensive treat- 
ment oi the region from this molecular line perspective. 
They estimate 4xl0 6 M of molecular gas, at the large 
end of GMC masses in the Milky Way. Thus, despite 
the prior star formation, there remains a tremendous 
reservoir of gas pregnant with star formation, offering 
the oppo rtunity to study ma ny evolutionary stages. To 



Motte et al.| ( |2007[ ) have carried out a contin- 
ey at 1.2 mm with MAMBO further revealing 



this end, 
uum survey 

through dust emission the complex and hierarchical mor- 
phology in the region. They find that Cyg X hosts about 
40 massive pro tostars destined to be OB stars. As dis- 
cussed in § |5.4[ the more massive protostars are already 
forming (ultr a) compact H II regions, famous e xamples 
being DR2 1 QDownes fc Rinehart|1966| ) , W75N (|Wester-| 
hout 1958 ), and AFGL2591, the strongest submillimeter_ 



CO emission integrated over the line profile (W) is of- 
ten taken as a surrogate of the column density of molec- 
ular hydrogen and so should be correlated with BLAST 
emission. We note the good correspondence and explore 
this in the examples below. Radiation from massive stars 
can of course in time destroy CO and the detailed corre- 
spondence with dust emission. 

Neither dust continuum emission nor W is sufficient 
for describing the complete physical environment and ge- 
ometry of the star formation region. Additional insight 
can be gained by exam ination of the CO velocity cubes 
( Schneider et al.||2006| ). Where there is CO coverage in 
our BLAST map, most of the identified BLAST compact 
sources are correlated with CO emission features (within 
the CO clumps mentioned above), thus associating the 
dust emission with gas at a certain velocity. In § |3.6[ 
we provide a few examples using their 13 CO J = 2 -)> 1 
cube to emphasize that objects seen closely together in 
projection in a dust emission image can have quite differ- 
ent velocities. In principle the velocity associations could 
be use d to sort the sources w ith respect to distance. In 
Cyg X |Schneider et al. (2006) have argued that the main 
velocity systems are all at roughly the same distance, 
1.7 kpc. Nevertheless, the association with different CO 
clumps having differing velocities does indicate that the 
BLAST compact sources are in distinguishable environ- 
ments, and there certainly are some sources behind the 
main Cyg OB2 complex (§ |4.2|). 



3.6. Selected Regions 

In this section we show the BLAST compact sources 
in the context of the diffuse submillimeter emission and 
other tracers discussed above. They are marked in the 
figures with a 1.3 7 square, characteristic of their apparent 
size in the deconvolved image. In each example, starting 
with Figure [6j the left panel shows the BLAST 250 /im 
image. On this are contours of the 13 CO J = 2 ^ 1 
emission integrated over velocity ranges whic h highlight 



compact sources within the BLAST coverage. Table [3] 
lis ts the corresponden ces between the cores and clumps 
of |Motte et"alf ( |2007D and our BLAST sources. 



the CO clum ps discussed in Appendix C of Sch neider 
|et al.| ( [2006 [ ). Good correlations between BLAST emis- 
sion and CO are revealed. The coordinate system chosen 
for these enlarged views is Equatorial to facilitate com- 
p arison with the cutout r egions in the annotated figures 
in |Schneider et ah (2006). Star clusters, local s ource s of 
power, ionization, and pressure, are also noted (§ 3.1). In 
the right panel is the MIPS 24 /im image with contours 
of the 21-cm radio emission from the CGPS map shown 
in Figure [2] Motte et al. (2007) show pairs consisting 
of their 1.2 mm MAMBO image and the corresponding 
MSX 8.3 jam image, but with no contours of CO or ra- 
dio emission. The MIPS image here is more sensitive to 
point-like sources, many coinciding with BLAST sources. 

The regions are presented in order of the number of 
the DR H II region in the field, whose nominal position 
is marked with an arrow or filled triangle in the right 
hand panel. To locate the regions in the large overview 
image in Galactic coordinates (Fig. [2]), see the DR num- 
bers marked there. Details of the nature of the BLAST 
compact sources and their evolution are deferred to §§[4| 
and[U 

DR7. The first region highlighted is a field including 
DR7, as shown in Figure|6] Immediately outside the DR7 
H II region itself is a rim of BLAST emission containing 
several compact sources. The typical rim shaped geome- 



10 



Roy et al. 




40:30 b 



20:32 



20:28 



20:31 20:30 20:29 20:28 20:32 20:31 20:30 20:29 

RA (J2000) RA (J2000) 

Fig. 6. — Left: BLAST 250 /im image inc luding DR7 over l aid w ith contours of 13 CO J = 2 — > 1 emission integrated over velocity ranges 
chosen to highlight the clumps discussed by |Schneider et aL] ( |2006| ), Fig. B.4. Solid, dash, and dot contour lines c orrespond to the ranges 3 
to 9, —7 to —1, and —13.5 to —8 kms -1 , respectively. Circle represents position and extent R50 of star cluster ( |Le Duigou fc Knodl seder 
|2002l. Right: 24 /im MIPS map of the same region, with contours of the 21-cm radio continuum emission from the COPS (see Fig. 
BLAST compact sources are marked in each image by square boxes of 1.3' width. Note the MIPS counterparts to many BLAST sources. 



try of the DR7 region has been formed by the star clus- 
ter (LK09) now residing in the cavity. Interestingly, the 
DR7 H II region along with its cluster are not associated 
with the contours of CO shown. This H II r egion has 
a rec ombination line velocity of —40 kms -1 (Lockman 
| 1989[ ) and CO velocity —50 kms -1 beyond the velocity 



plex near kms -1 is be ing in fluenced by Cyg OB2 off to 
the upper left (see also 
IRDC ridge containin; 
As shown in Figure [3 



5.1.2). This range also traces the 
U7^34+0.33 and G79.27 +0.38 . 

detail in 8 15.6.2 



coverage i n|Schneider et al. ( 2006| ) , placing it well behind 
the Cyg OB2 complex, at 3.6 to V.5 kpc, possibly in the 
Perseus arm. Another indicator of the larger distance 
is the lack of signature in the extinction map, which is 
based on colors of detectable 2MASS stars. 

There are, however, several CO clumps that do coin- 
cide with BLAST emission and compact sources. The 
strong BLAST emission (with compact sources with 
MIPS counterparts) to the south coincides with clump 
5. Not all clumps (which have different W) are easily 
seen in the BLAST emission (e.g., clumps 7, 8, and 9). 
The bright BLAST source at the top, in clump 1, has 
a compact H II re gion and a cometary tail linking it to 
Cyg OB2 (§[5T2|. 



The emission at the upper left, associated with clump 
2, is perhaps the most interesting, having a chaotic mor- 
phology in both MIPS and BLAST images (see also 
Fig. m). In the ridge we find three BLAST sources. Al- 
though in projection the center of Cyg OB2 is close by, 
about 15' to the east (left), there is no sign of interaction; 
furthermore, there is a parallel ridge of H II emission on 
the right, suggesting illumination from that side. It co- 
incides with G79. 957+0. 866 for which the re combination 
line velocity is —14.8 kms -1 (Lockman 19891), close to the 
CO velocity of —11 kms -1 . In the entire map, this CO 
cloud is unique in this velocity range < — 10 kms -1 and 
not obviously connected to the other Cyg X molecular 
complexes. This then seems to be an object so mewha t 
more distant, beyond the influence of Cyg OB2 (§ |5.1.2 ). 

DR15. Figure [7] shows the fascinating region contain- 
ing DR15 an d the above-mentioned PAH arc. Schneider 
et al.l (|2006|) argue that the DR15 molecular cloud com- 



and in more detail m 
BLAST sees this cold ridge in emission and resolves sev- 
eral compact sources; in the MIPS image, some embed- 
ded protostars shine faintly through. Toward the center 
of this field is the protostar IRAS 20293+3952, part of 
a different CO system at positive velocity. The pair of 
BLAST sources to the upper right are associated with 
streamers in the MIPS image which appear to be "blown" 
from the right. 

DR17N. The lower portion of Figure [8] shows the north- 
e rn pa rt of the extended DR17 region (see also Fig.[2]and 
§ |5.4[) . A major arc-shaped extended structure crosses 
BLA&T and MIPS images, and is seen in 8 /im PAH 
emission as well, confining the H II emission influenced 
by the OB clusters inside DR17. CO in the higher ve- 
locity range 6 to 14 kms -1 best traces the arc. The 
lower velocity system projects across this, and has its 
own BLAST emissi on and compact so u rces 



DR17-Pillars. Schneider et al. (2006) identified 



"molecular pillars" (their Fig. B.l) in the higher velocity 
range whose orientation points to the influence of the OB 
stars in DR17. BLAST finds compact sources associated 
with these pillars as shown in Figure [9j Compared to 
other dust in DR17, the above-mentioned arc and these 
molecular pillars are cool (Fig. [3]). There is another H II 
region to the left in the image, with BLAST emission 
and compact sources alon g its interface as in DR7. This 
is the "Diamond Ring" ( Marston et al.|2004 ) at an inter- 
mediate velocity (8 kms -1 ), as described further in the 
DR21 discussion below. 

DR20. Figure [10] shows DR20, which is at the end of 
a prominent ridge seen in the low- velo city CO (clump 1) 
and BLAST emission (see Fig. 13 in § 4.1 below). There 
are several BLAST compact sources here and in the com- 
plex to the west (clump 2, DR20W). The source in the 
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Fig. 9. Same as Fig.^but for the DR17-Pillars region to its south ([Schneider et al.|2 006). Solid and dashed CO contours are for ranges 
1 to 6 and 14 to 20 kms , as used in Fig. B.l of Schneider et al. (2 006^ The long molecular pillar on the right is that labelled in Fig. 3 
A non-thermal extragalactic radio source 18P 69 ( Wendker et al.[199l|> also prominent in Fig. [2] is marked. 
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RA (J2000) RA (J2000) 

Fig. 10. — Same as F ig. [6] but for DR20. Solid and dashed CO contours for ranges —10 to 1 and 6 to 14 kms -1 , as used in Fig. B.2 of 
Schneider et al.| ( |2006| >. They called the molecular clump to the rig ht DR 20W, in the same velocity system as DR20, and that containing 
m a different velocity system of a prominent elephant trunk (§ |5.1.1[ ) DR20NW. 
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middle of this image is associated with higher velocity 
gas (clump 3, DR20NW), as is more diffuse emission and 
other compa ct sources along the pillar toward the upper 
left (§ |5.1.1|) , possibly connected to the velocity system 
of the DK17 molecular pillars and arc. 

DR21 and W75N. The DR21 and W75N regions shown 
in Figure [TT] is obviously very active. The ridge of 
BLAST emission including DR21 and DR21(OH) points 
deceptively to W75N. However, the peak CO emission 
on the ridge (clumps 3 and 4) and clump 8 containing 
the BLAST sources to the north west of W75N are at 
—3 kms -1 (the low velocity system), whereas clump 10 
of W75N is in a quite different higher velocity compo- 
nent at +10 kms -1 . The low velocity system is fairly 
widespread, including a complex to the south west (over- 
lapping the lower edge of the eastern-most DR17 molec- 
ular pillar) and extending down to DR23 and DR22 
(Fig. 12). The mass column density traced by BLAST 



correlates well with the CO emission. A fascinating fea- 
ture is the extension of the DR21 ridge to the south, 
which shows up as a very narrow dark lane in the 24 /im 
image. Th is line ar feature is at —4.5 kms -1 (see also 



Fig. 28 in § 5.6.2 below). 



Thenigher velocity gas also contains prominent sources 
other than W75N, inclu ding the "Diamond" (clump 5) 
of the "Diamond Ring" (jMarston et al.|2004| ). The more 
extended H II region along the "Ring" appears to be 
defined by both velocity systems. 

DR22 and DR23. The CO cloud structures in the 
DR22 and DR23 regions shown in Figure [12] correlate 
well with the BLAST dust emission. Immediately to the 
south west of DR23 is a dim bay in the MIPS image, 
where the main CO cloud (clump 1) of the low- velocity 
system confines the H II region. Along this interface 
are three BLAST compact sources, including C132, 
perhaps triggered by the earlier generation of stars 
in the DR23 cluster DB17. The two DR regions are 
interconnected by a CO filament (several clumps) in 
the lower velocity range. Along this filament there are 
BLAST compact sources and dark clouds, one (clump 3) 
wi th a p rotostellar nursery visible with MIPS and IRAC 



(§ 15.6.21). For DR22, the complementar y det ail provided 
TAC 



by IRAC emission (see Fig. 23 



5.2) confirms 



classical PDR geometry, with the lower left being most 
edge on. The location of the BLAST sources, not just 
the compact H II region but all around the periphery, is 
suggestive of triggering. 



4. PROPERTIES OF COMPACT SOURCES 

4.1. Submillimeter-MIR SEDs 

The SED for cold dust emission at temperature ~ 13 K 
and emissivity index f3 = 1.5 peaks at 250 /im, and so 
the combination of the three BLAST filters alone is well 
suited for determining the dust temperature. Neverthe- 
less, there is a range of temperature among sources (§ 4.3 ) 
and it is always preferable to have a broad wavelength 
coverage spanning both sides of the peak. Fortunately, 
Cyg X is well covered by both blind and more-targeted 
surveys. An example of the mult i- wavelength coverage is 
given in Figure [13 Clearly, resolution changes from map 



to map. Also, for the central source here, the centroid 
changes at sh orter wavelengths (and in the radio). As 
described in § |2.5[ we have extracted flux densities con- 
sistently within the same size beam, to characterize the 
same physical structure. 

We fit an idealized single-temperature SED expressed 
by equation ^ to data at all available wavelengths to 
determine tempe rature, mass, and to tal l uminosity of 



2008 



each source (see |Chapin et al.| 
|2008| for details). We adopt th 
ues, namely = 10 cm 2 g _1 , r 
and D = 1.7 kpc except for a few distant sources ( 



and Truch et al. 
e same parameter val- 
= 0.01, and f3 = 1.5, 

Single-temperature SEDs based on this value of f3 fit the 
data of typical sources well. One of the main systematic 
uncertainties is the value of K$r which is probably not 
known within a factor of two. 

We treated ancillary photometric data for wavelengths 
less than 100 /im (see AppendixfA| in the specific context 
of each individual SED, including them as upper lim- 
its where available/necessary. Upper limits come with 
an additional penalty functi on in the y 2 minimization 



through a survival equation ( |Chapin et al. 2008 Truch 



et al 



2008). 



BLAST filters have large spectral widths about the 
central wavelength. Color correction to monochromatic 
flux densit y is carried out as p art of the SED fitting de- 
scribed by Truch et al. (2008). The factors are close to 
unity and fairly consistent from source to source for the 
range of temperatures encountered here: 1.04, 1.07, and 
0.99 for 250, 350, and 500 /im, respectively. 

An example of a multi-wavelength SED fit is shown 
for C169 in Figure [l4j this is one of the most luminous 
sources, W75N, andone of the hottest, with a best fit 
temperature for mm-FIR data of 36 K. In contrast to 
this, the SED of one of the colder sources is shown in 
Figure [15] this is C116 with a temperature of 17 K. 

In Table [4] we record T, M, E, and L and their uncer- 
tainties for 170 individual sources, there being no entries 
for the 14 faint sources either near the map edge or with 
unreliable 500 /im flux densities. The uncertainties in T, 
M, and L, and the corresponding 68% confidence enve- 
lope of possible modified blackbodies, we re obtained by 
the Monte Carlo technique described by |Chapin et al. 
( [2008| . 

4.2. Sources Behind Cyg OB2 

There are eight sources at distances larger than the 
main Cyg OB2 complex, C27, C30, C33, and C35 in DR7 
(§|3^), C59, C60, and C61 in G79.957 + 0.866 (§§ HT6| 
anc 



.1.2), and C71 in ECX6-27 (§§ 2.1 and 3.1). 



calculating L and M we assigned these rough distances 
of 6.8, 3.4, and 8.5 kpc, respectively, making L and M 
16, 4, and 25 times larger than if the Cyg OB2 distance 
of 1.7 kpc were adopted. This makes C30 the most lu- 
minous source and C71 the most massive. Both are of 
course integrated over much larger spatial scales than 
for the sources at 1.7 kpc. Note that L/M and T are 
preserved. 

4.3. Temperature 

Temperature is obtained directly as a free parameter 
of the SED fit. Figure [16] shows the temperature his- 
togram of BLAST sources detected in the Cyg X field. 
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Fig. 11. — Same as Fig, jglbut for D R21 and W75N regions. Solid and dashed CO contours are for ranges —10 to 1 and 6 to 14 kms 
as used in Fig. B.3 of Schneider et al.|fl2006) . 



The superimposed h istogram is for the Vulpecula region 
( Chapin et al.||2008 ). The median of the distribution in 
Cyg X is about 22 K which is comparable to the median 
for Vulpecula at 21 K. Recall that these are for j3 = 1.5. 
A somewhat higher (3 combined with a temp erature dis- 
tribution might be a more realistic model ( Netterfield 
et al. 11 2009]) , though in the simple model here there is a 
T~F< 



massive stars, as o ppo sed to fragmentation into much 



lower masses (see § 6.3 for further discussion). 



The surface column density for fragments of radius R 
is given by 



2.7 x 10" 



single T. For /? 
10% lower. 



2, the temperatures would be about 



M 



100 Mrr 



R 



0.5 pc 



g cm 



(3) 



4.4. Mass 

Source mass M is obtained from the best fit amplitude 
and temperature of the SED and equation (J2|, which 
shows how any uncertainties in ft, r, and D are propa- 
gated as systematic errors. Recall that our sources typi- 
cally have an apparent FWHM of 1', which corresponds 
to about 0.5 pc. This is more characteristic of a "clump" 
mass, as distinguished fro m the h i gh res olution measure- 



ments of "cores" by |Motte et al i 

values adopted for kq and fj do provide cons istency with 
the 1.2 mm opac ity of 1 cm 2 g -1 adopted by Motte et al 



We adopt the deconvolved FWHM of the brightness pro- 
file as the radius. Surface densities of the cores (and 
clumps) of |Motte et al.|(|2007[ ) are listed in Table|3] (Note 
that the assumption of a temperature of 20 K for the 
cores means that the surface density reported might be 
an upper limit.) As discussed below, the surface density 
of BLAST clumps is somewhat lower, in part because 
of the low angular resolution which limits the ability to 

discern sub-structure at small R. 

Following equation (2) of |Motte et al] Q2007|), the 



( 2007J) . (Note that the volume- averaged molecular hydrogen density is 



( [2007] ) , following |Ossenkopf fc Henning| Q1994p . H owever, 
without the multi- wavelength coverage, they were forced 
to adopt a temperature.) From these derived masses, 
ranging from 10 to 2000 M , the BLAST sources have 
the potential to form (many) massive stars under favor- 
able physical conditions. For example, in the context 
of the turbulent core model |Krumholz fc McKee ( 2008 ) 
predicted a critical surface density tor the formation of 



(n H2 ) = 3.3 x 10 3 M/(100 M ) cm" 



(4) 



Because of the BLAST beam size, these clump densities 
are more than an order of magnitude lower th an for the 
embedded cores found by Mo tte et al. l (j2007|). The c or- 
responding free-fall time scale (Stahler fc Palla 2005]) is 
also longer, 



t s = 5.6 x 10 5 [M/(100 M 



-1/2 



yr. 



(5) 
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Fig. 12. — Same as Fig. [n] but for DR22 and DR23 regions to its south. 



4.5. Luminosity 

The best fit single-temperature SED is a good and inte- 
grable interpolating function from which the far-infrared 
luminosity L is derived analytically. The frequencies at 
which the most important contributions arise can be vi- 
sualized most readily when the SED is plotted in vF v 
form like Figure [T4| 

Often the short- wavelength mid-infrared data lie in ex- 
cess of this simple SED and to allow for that we sim- 
ply integrate the piecewise curve connecting that short- 
wavelength data where available. This is normally a 
small contribution to the total bolometric luminosity 
L bol (Table [§. For example, for W75N (C169, Fig.[l4) 
L rises from 6.2 to 6.6 x 10 4 L . In the context of posi- 
tion in the logarithmic L-M plot below, this is not a big 
correction. 

The luminosity range in Cyg X is large, extending over 
three orders of magnitude from 7x 10 4 L down to 40 L , 
beyond which we lose many sources in the (cirrus) noise 
(note that the a pproximate 250 /im completeness line is 
about 30 Jy; see Chapin et al.||2008 ). 



4.6. The L-M Diagram 
The L — M diagram can be exp loited to assess evolu- 



tionary stages (Molinari et al. 2008). Figure 17 shows our 
results in the L — M plane lor those sources at the dis- 
tance of Cyg OB2; the results for the few distant sources 
(§|4.2|) are in Figure [l8| (note the scales to larger masses 



and luminosities). 

Lines of constant T are diagonal lines (constant L/M) 
in this L-M diagram: a mass M radiates a predictable 
L, depending on n^r and /3. To be consistent with the 
analytical loci, we plot L from the SED fit rather than the 
only slightly larger Lboi- "Orthogonal" to these diagonal 
lines are loci of constant 250 jam flux density. Note that 
there are relatively fewer sources between the 30 and 
10 Jy loci because of the growing effect of cirrus noise. 

In this diagram the "error ellipse" determined by the 
Monte Carlo technique is elongated along the locus for 
the source flux density and the extent is well described 
by AT/T. We examined the histogram of AT/T and 
rejected outliers > 0.19; these were 18 sources with 
poorly constrained SEDs, usually lacking in definitive 
data near 100 /xm combined with their poor image qual- 
ity at 500 jam. 

The dashed and dot-dash thick lines, roughl y lines of 
constant T, are the loci obtained empirically by Molinari 
|et al.| ( |2008[ ) for sources thought to be in the accretion 
stage and the later nuclear burning stage (when envelope 
dispersal begins), respectively. For single low mass stars, 
these would correspond to the spectroscopically defined 
"class 0" and "class 1," respectively. With BLAST, we 
detect sources of the size of "clumps" and furthermore 
even smaller angular-size high-mass "cores" are capable 
of forming multiple stars. This warns against a simple 
interpretation of this diagram based on single-star evo- 
lutionary tracks, although it is possible that once high 
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BLAST 500 fJ,m BLAST 350 /*m BLAST 250 /xm 



IRAS 100 jtxm MIPS 70 /xm IRAS 60 fj.m 



MIPS 24 /xm MSX 8 /xm CGPS 21 cm 



Fig. 13. — Thumbnails zooming in on a 14' by 18' section of 
the BL AST survey area containing DR20 (Galactic coordinates; 
cf. Fig. \10\ at multiple wavelengths available for photometry. Rel- 
ative appearance of structures changes because of different dust 
temperatures. The lowest in the vertical chain of three sources, 
seen only at wavelengths 500 jum and longer and so not catalogued 
as a BLAST source, must be quite cold. At IRAS 60 /im (IGA, 
after HIRES processing), sources appear elliptical across the scan 
direction. At IRAS 100 /im, the emission of crowded sources is 
often blended. 

mass stars form, the most massive will dominate the lu- 
minosity and ionization (this will depend on the IMF, 
star formation efficiency, and also small-number statis- 
tics). 

We think that the most illuminating way to think of 
pre-stellar evolution in this diagram is in terms of the 
energy source for the clump, which determines the ap- 
propriate equilibrium temperature T for the approximate 
SED. In the very earliest stages being sought in submil- 
limeter surveys, the energy source for the clumps is pre- 
dominantly external, namely the impinging interstellar 
radiation field. In Cyg X the radiation field is in princi- 
ple higher than in the local interstellar medium because 
of all the massive stars that have already formed, but 
this radiation is attenuated by the dusty molecular ma- 
terial in which the pre-stellar clumps are embedded. A 
useful point of reference is the equilibrium temperature 
corr espo nding to L = M, which for the adopted n^r and 
P (§[4TT1) is 16.2 K. 
Massive sources located below this L = M locus (see 



5.6) still can have substantial L by virtue of a lot of 



mass, but cannot have any predominant internal source 
of energy (either accretion or nuclear), and so could be 
called "starless." This term is probably best avoided for 
these clumps, within which there might actually be a few 
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Fig. 14. SED of C169 (W75N). To show intuitively where 
the most important contributions to the bolometric luminosity 
arise, this particular log-log plot uses vF u , obtained by multiply- 
ing flux densities by v/v2b0 with V250 corresponding to 250 fim. 
The central solid curve shows the best-fit modified blackbody us- 
ing data at A > 60 jum with f3 = 1.5. The MSX data shown 
by the circles ar e no t used in the fit but are important in con- 
straining Lbol (§ |4.5| ). The bracketing curves represent the 68% 
confidence envelope of modified blackbody models obtained from 
Monte Carlo simulations. Best fit parameters are T = 36. 2 ± 3.6 K, 
M = (7.6 ± 1.2) x 10 2 M , and L = (6.2 ± 2.2) x 10 4 L . 
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Fig. 15. — Like Fig.[l4]but for one of the colder sources in Cyg X, 
C116. At 70 /im we plot the 1-cr (arrow) and 3-cr upper limits 
(Appendix [Al; upper limits c onstrain the SED through a penalty 
function (CiTap in et al.|2008| >. Best fit parameters are T = 17.1 ± 
0.9 K, M = (1.7 ± 0.2) x 10 2 M , and L = (2.2 ± 0.4) x 10 2 
Lq. The 24 /im MIPS image in the upper panel of Fig. |26| reveals 
a deeply embedded stellar nursery but it is not yet very luminous. 
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Fig. 16. — Light grey histogram shows the distribution of source 
temperatures for the Cyg X field. Dark over-plotted histogram is 
for the BLAST05 Vulpecula field jChapin et al.|2008| >. 

low mass stars already formin g an d detectable by sensi 



tive telescopes like Spitzer (§ |6.2| ). The key considera- 
tion is what is the dominant source of energy determin- 
ing T. As mentioned in the introduction, these clumps 
could be said to be in stage E ("E" for "external" or 
"earliest" ) , which seems to us a better terminology than 
"starless," or "class —1" for single-star pre-stellar cores. 
Figure [17| shows that there are many clumps in this stage. 
We can calculate the Bonner-Ebert mass, above which 
(within the assumptions) a clump is gravitation ally un- 
stable (e.g., |StahlCT^P^[200^: 

M BE = 1 [T g /(10 K)] 3 / 2 [n H2 /(10 4 cm" 3 )]- 1 / 2 M . (6) 

Thus, assuming that T g < T, the low-luminosity, low- 
temperature BLAST clumps are unstable. As discussed 
below, in the simplest theory they would evolve into 
the higher-luminosity, higher-temperature clumps seen 
higher in the diagram at the same M. 

As gravitational collapse progresses, more and more 
mass accretes into protostellar cores. The actual lumi- 
nosity in this accretion-powered stage depends on the 
accretion rate and the potential well. In recognition of 
the underlying energetics, these clumps could be said to 
be in stage A ("A" for accretion-powered). Consistent 
with this interpretation, sources observed to be in this 
part of the L-M diagram have a characteristic signature 
of active accretion (§ 5.5). 

Ultimately, nuclear Fusion becomes the dominant 
source of power. Unlike for low mass stars, massive 
pre-main sequence stars probably continue accreting af- 



FlG. 17. — Distribution of BLAST sources in the Cyg OB2 com- 
plex in the L-M plane. Those to the right of the steep dash-dot 
(gray) curve are above the Bonner-Ebert mass. Solid diagonal lines 
are loci of constant T or L/M. The dotted line (red) is L = M. 
Dot-dash curves "orthogonal" to these are for constant 250 /im flux 
density. Dash (green) and dot-dash (blue) lines denote the location 
of sources powered by accretion and nuclear burning , respe ctively 
as derived empirically in Fig. 9 of |Molinari et al.| ( |2008] >. Cool 
sources discovered by BLAST with low L/M still are externally 
heated (stage E). These appear to be gravitationally bound, but 
have no significant internal power from star formation yet. The ver- 
tical arrow indicates the direction of evolution in this diagram as 
protostar formation takes hold within a clump. The horizontal ar- 
row indicates the direction of evolution as the embedding material 
is dispersed by the formed stars/cluster, assuming the surrounding 
dust still reprocesses most of the internally-generated L. Other- 
wise (as appears to be the case), the re-radiated L, measured here, 
is less. B lack filled circles a re BLAST sources corresponding to 
clumps of Motte et al. (2007) (Table [3}. Crosses are sources with 
morphological evidence ol mass-stripping from radiative interac- 
tion with Cyg OB2; some at lower luminosity exhibit the effects of 
external ionization. 



ter first beginning nuclear burning , increasing their mass 
further ( jZinnecker fc Yorke||2QQ7 ). When this accretion 
ends, an individual star is at its final position on the 
ZAMS. It will still be embedded, and so optically ob- 
scured, but its vast power reprocessed will produce a rel- 
atively warm far-infrared source. For sufficiently massive 
stars, significant ionization of the surroundings will en- 
sue, producing a hypercompact H II region. The re could 
be many within a sin gle clump, for example ( | Rivera- 1 
Ingraham et al. 2010). These clumps would appear m 
the L — M diag ram near the empirical nuclear burning 



locus (see § 5.4). 



As the massive stars clear their local environment 
through the expansion of the H II region and via ra- 
diation pressure on dust (successively ultracompact then 
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Fig. 18. — Li ke Fig. [TtJ but for sources more distant than 
Cyg OB2 (§ |4.2| ). Triangles, circles, and square represent sources 
about 3.4, 6.8, and 8.5 kpc away, respectively. These include the 
most luminous and most massive sources in the survey. 

compact H II regions), the individual objects will become 
more extended far-infrared sources. The strong short- 
wavelength stellar radiation field will enhance both ther- 
mal emission and non-equilibrium emission from PAHs 
in nearby photodissociation region material, making the 
SED of the integrated re-emission of the clump broader 
and more complex. Evolution in the L-M diagram will 
initially be at constant L with decreasing M; ultimately, 
the optical depth and/or covering factor would decrease 
so that the reprocessed L will fall. 

A stellar cluster with an extended H II region would 
emerge. Even though such clusters are found in Cyg X, 
they are not detected as compact BLAST sources and 
so do not appear in the L-M diagram; only objects like 
their precursors would. 

4.7. Census 

It is interesting to ask how many objects are in each 
evolutionary stage. To address this we have concentrated 
on the mass range between 50 and 500 M where there is 
an abundance of sources (100 omitting the sources more 
distant than Cyg OB 2) and the lower luminosity end is 
minimally affected by completeness, at least down to L = 
M. Figure [l9| presents a histogram showing the relative 
populations in terms of L/M, which we have argued is at 
least qualitatively related to the successive stages. Note 
that this figure looks quite similar to Figure [l6j there 
being a non-linear mapping of L/M into T. 

Figure [19| shows that there are many sources in stage 
E (L/M < 1 L /M ), even though it is clear that 
there are selection effects (basically the combination of 
low flux density and cirrus noise) beginning to limit their 
detect ability. Most of the sources are in the accretion- 
dominated stage A (L/M up to about 30 or T = 30 K). 
The relatively fewer hotter sources, where nuclear burn- 
ing is taking over, presumably reflects the shorter lifetime 



c 



CD 
_Q 



30 



20 



10 



Stage E 



Stage A 



-2-10123 
Log (L/M)/(L /M ) 

Fig. 19. — Histogram of L/M ratio of the Cyg X compact sources 
in the mass range 50 - 500 Mq . The dotted vert ical line separates 
stage E sources below the locus L = M in Fig.j^] Most sources are 
in the accretion-powered stage A, up to the dashed line (L/M « 
30 L0 /Mq), beyond which nuclear burning becomes predominant. 

of this stage, which could be characterized by more rapid 
final collapse and then envelope dispersal. 

5. EVOLUTIONARY STAGES 

Although our BLAST survey is "blind and unbiased", 
in surveying the entire region rath er than selecting sub- 
regions with, say, high extinction ( jSchneider et al.||2QQ6| 
Motte et al. 2007), it is not unbiased in another sense. 
This region is obviously well known for its GMC and 
having formed the Cyg OB2 association, which has dra- 
matically influenced the surrounding molecular material 



(§5.1). In such a region one expects there to be both 
triggered and/or sequential star formation, supplement- 
ing spontaneous star formation. The conditions could be 
quite different now than what preceded the formation of 
Cyg OB2. We find no evidence for the precursor of an- 
other such massive compact association. Nevertheless, 
within this reservoir several smaller, but sti ll notable, 
embedded clusters have been identified (§ |3.1|) . The OB 
stars in the more evolved ones hav e produced extended 



H II regions like DR17 and DR22 flDownes fc Rinehart 
1966). 

13 CO 



As discu ssed in S |3.5l in thei r 13 CO J = 2 H> 1 
data cubes, [Schneider et al.| (|2006) have identified many 
clumps (detailed in their Appendix C). Where the cov- 
erage overlaps, each of the BLAST compact sources can 
be linked through t he m orphology to one of these clumps 
(see examples in § 3.6). In Figure 20 we c omp are the 
sum of the masses of the BLAST sources (§ |4.4[) within 
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compact precursors to such clust ers m i ght be detectabl e 
in the Planck Cold Core survey (§ [Z2||Juvela et al.|2010[ ). 
It will of course depend on how long-lived the precursor 



5.0x10 3 1.0x10 4 1.5x10 4 2.0x10 4 2.5x10 4 3.0x10^ 
CO Clump Mass [M ] 

Fig. 20. — Sum of the m asses of the BLAST sources within CO 
clumps of Schn eider et al.| ( |2006| > versus the CO-estimated clump 
mass. Linear correlation has a slope of 0.05. An outlier at 3000 
Mq along the y direction, pointed to by the arrow, is not included 
in the fit. 

each clump to the tabulated CO-estimated clump mass. 
These are fairly well correlated, and we find typically 
that only a few percent of the material over these ex- 
tended regions is in the form of compact sources, in some 
way related to the (potential) star formation efficiency. 
This plus the timescales for clump evolution indicate that 
there is still an interesting future for star formation in 
this GMC. 

From this perspective on the complex range of star for- 
mation going on in Cyg X, our goal is to work toward 
submillimeter evidence for the earliest stages of massive 
star formation. The somewhat later stages have been de- 
tected by a variety of earlier observations, and BLAST 
sees these too. The many clusters in Cyg X are agents 
producing complex morphologies in the medium, through 
expanding H II regions and strong UV radiation fields, 
and they might play a pivotal role in inducing further 
star formation in the molecular clouds. We work in re- 
verse chronological order, finding what the outcome of a 
particular stage looks like in multi-wavelength data, and 
then asking what this tells us to look for in the stages 
that came before. 



5.1. The Influence of the Massive Cluster Cyg OB 2 
Comparing the estim ated mass of Cyg OB2 ( 
4 to 10 x 10 4 M dKnodlseder]|2000|) to the r emaining 



3.1), 



molecular reserve in Cyg X, 4 x WM Q (§|3.5l Schneider 
et al.|2QQ6| ), indicates that formation of this Un2 associa- 
tion was the major event for this GMC, perhaps not to be 
repeated. Both of the central clusters, BBD1 and BBD2, 
are more massive than the Trapezium. Again, will there 
be more elsewhere in Cyg X? There is no BLAST clump 
that could be a precursor to such a prominent cluster 
and even adding the mass of an embedding CO clump 
(the most massive of which is 2 x 10 4 Mq) would fall 
short. Nevertheless, in unevolved GMCs, massive and 



depend on now long-lived tne precursor 
stage is. 

Massive young stars, through ionization and radia- 
tion pressure, dramatically influence their e nvironment 
Cyg OB 2 h as an age of about 2.5 x 10 6 yr ( |Negueruela 



et al.|[2QQ8| ), and in that time has created in the ISM 
a lower density region of radius 1° (30 pc), ionized in 
the interior (see Fig. [2]) and surrounded by distinctive 
molecular cloud complexes. The connectivity of these 
clouds in the CO data cube and the signs of interaction 
through UV rad iation (MSX Band A data) led Schnei- 
der et al. (2006) to conclude that they were all located 
at the same distance. While obviously not yet disrupted, 
the clouds have been reorganized and so Cyg OB2 might 
have induced some of the further star formation that is 
on-going. 

Above Cyg OB2 in Figure [2] is a relative void (nev- 
ertheless, it is projected on an emission and extinc- 
tion plateau), with no high-contrast structures of sig- 
nificant column density in any tracer (molecules, dust, 
ionized gas). The right hand edge is defined by CXR9 
whereas the left, though well-defined and containing 
some BLAST emission too, is unnamed. This void could 
be a "chimney" blown out by Cyg OB2, channeling high 
pressure gas out of the Galactic disk. 

5.1.1. Pillars 

Radiative forces sculpt pillars, where there are pre- 
existing dense molecular structures. The range of radia- 
tive influence of Cyg OB2 extends to a radius of at least 
2° (60 pc), as illustrated by the prominent pillars near 
G79.4-1.1 and G81.4+2.2 in the Spitzer IRAC and 24 /im 
MIPS images. Figure [21] shows an example of a long pil- 
lar or elephant trunk pointing toward the Cyg OB2 core. 
At the end of the trunk is a luminous BLAS T clump 
C85. In CO em ission, it is called DR20NW by |Schnei- 
| der e t al" ( 2006[ ), with CO velocity +12 kms -1 (see also 
L'he clump is the obstruction defining the trunk 



structure. Star formation is occurring in this clump, sug- 
gestive of triggering by Cyg OB2. This geometry is seen 
elsewhere; for example, C75 on the left edge of the chim- 
ney is also at the head of a pillar facing Cyg OB2 (see 
also the cometary tails below). 

Such a peninsular structure is continually eroded, 
which would produce eventually an isolated clump with 
a cometary tail, but this particular trunk appears to 
be attached to the molecular cloud, curving round and 
broadening out to include C91, C98, and C100 along the 
interface. These structures also evolve due to instabili- 
ties at the interface of the H II region and the neutral 
molecular cloud, and evidently on-going star formation 
is commonly induced. 

Note that DR20 (C88, C87, C90) is at a projected dis- 
tance of 6' from this trunk, and yet does not display the 
same hallmark interactions. However, this complex is at 
a very different velocity —3 kms -1 , mo re closely associ- 



ated with DR21. |Schneider et"aT| ( 2006[ ) note that clump 
DR20W, at the same velocity, does have an elongated 



shape pointing to Cyg OB2 (see Fig. 10). BLAST C80 



is at the brighter head of its CO emission and C83 in 
the dimmer tail. Inspection of the Spitzer images shows 
signs of interaction with nearby stars in Cyg OB2, but 
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Fig. 21. — IRAC band 4 image showing a prominent pillar or 
"elephant trunk" structure pointing to Cyg OB2. Squares are 
BLAST sources and stars show the position of the OB stars. 

at a different position angle than the CO clump which is 
o riented more toward the core. 

|Schneider et al. (2006) also identify other molecular 
pillars in the CO channel maps, the most prominent 
of which is the long "DR17 pillar" at a velocity about 
— 13 kms -1 (Fig.JM. From its direction of elongation, it 
appears to be influenced by the cluster LK12 in DR17, 
but there is not perfect alignment (possibly there is some 
influence of LK14 to the north). This cluster does pro- 
duce a bright mid-infrared interface at BLAST CI 15, 
near the end of the pillar. The molecular pillar stretches 
from C115 through C114.C112, C113, and C108 to C109 
and possibly C106 (Fig. [9]). The CO channel maps hint 
at a connection between this DR17 pillar and the above- 
mentioned elephant trunk, but this pillar does not appear 
to be influenced by Cyg OB2. 



5.1.2. Cometary Tails 

Further striking evidence for interaction is provided 
by cometary tails in more detached structures, pointing 
away from the core of Cyg OB2. These too are seen most 
clearly in Spitzer IRAC and 24 /am MIPS images. 

DR18 is a large scale prototype, with a ionized edge 
facing the Cyg OB2 core, the BLAST clump C73 in the 
neutral region next to this, and an extended tail seen 
in both CO and dust emission (mid-infrared to submm). 
It is actually more complicated than this in detail, be- 
cause a loose aggregate of early B stars has formed (trig- 
gered?) near the leading edge adding additional ionizing 
power and sculpting of the immediate PDR environment 
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Fig. 22. — 24 /im MIPS image showing two BLAST sources, 
C38 and C17, with cometary tails. The spatial extents of the tails 
are roughly 1.5 and 3 pc, respectively. Contours of 21-cm radio 
emission from CGPS show ridge CXR11. 



flComeron fc Torra)[l999l ) . 

Another interesting example in the opposite direction 
is the cluster ECX6-21 (DB12) containing BLA ST C55 
in the head. C55 corresponds to OB2 globule 2 of[Schnei- 
der et al.| (2QQ6| with CO velocity -4.5 kms" 1 . About 
5' "downstream is a miniature version containing C53. 
These are located on the inner edge of the ionized ridge 
CXR9. 

Further afield from Cyg OB2 is DR15, with a cluster 
and two BLAST clumps, C32 and C28. Here the puta- 
tive tail (see Fig. [7]) appears to be part of an extended 
structure prominent in presumed PAH emission, which 
bends significantly toward lower latitude at a distance of 
17', near C12. There are several other BLAST clumps 
along this structure, C21 before the bend, and Cll, C9, 
C4, and C3 beyond. 

Figure [22] shows a field quite close by, spanning the 
ionized ricTge CXR11. There are two BLAST clumps 
C38 (on the left) and C17 with material ablated from 
the clump clearly being blown away. However, in the 
BLAST image these cometary tails are not detected be- 
cause of low column density. C38 is massive enough to 
have a detectable ionized edge. The C17 tail extends to 
an angular distance of 6' or about 3 pc in spatial scale. 
If, for illustration, the clump lifetime were 10 5 yr the re- 
quired material speed would be 30 kms -1 to reach such 
an extent. 

There are many other sources like this with obvious 
interactions. Their positions in the L-M diagram are 
highlighted in Figure [T7| For these exposed sources, it 



is clear that a sizable amount of mass is lost due to this 
erosion, and so the more evolved sources occupy a region 
towards the lower end of the mass axis in the L-M plot. 

The ridge of submillimeter emission containing sources 
C59, C61, and C60 provides an interesting counter- 
example. The superimposed CO ridg e is called "OB2 
globule 1" by |Schneider et al.| ([2006 ) who suggest an 
interaction due to its compactness and proximity to the 
core of Cyg OB2, in projection about 15' away. However, 
inspection of the detailed Spitzer images now available 
reveals no signs of interaction. Furthermore, the dis- 
placed ionized ridge parallel to the dust emission is on 
the side away from Cyg OB2 (see Fig. [6|. Together with 
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Fig. 23. IRAC band 4 (8.6 /xm) image of DR22. Ci rcle marks 
the position and siz e (R50) of the LK11 OB star cluster ( |Le Duigou] 
|fc Kno dlseder 2002). Overlaid contours are 21-cm radio continuum 
emission irom the CGPS. Squares indicate the BLAST sources, 
embedded in a rcs of submillimeter emission outside the ionization 
front (Fig. [12). 

the peculiar velocity, the evidence is that this object is 
more distant, beyond the influence of Cyg OB2. 

5.2. Clusters and Extended H II Regions 

Interactions of clusters with the ISM produces large 
scale complex structures and directly affect the next 
generation star formation efficiency in the region. The 



Cyg X region hosts numerous young open clusters (§ 3.1 ). 
In the BLAST region, there are five which have well- 
resolved H II regions with a fairly classical geometry. 
Expansion of the ionized gas, stellar winds, and radiation 
pressure have evacuated the core of the parent molecular 
cloud and produced a dense shell which appears in an 
arc-shaped geometry surrounding the cavity and cluster. 
Partial shells are suggestive of a blow-out or blister ge- 
ometry. The 21-cm radio emission traces the ionization 
front, and slightly outside of this is the FIR and submil- 
limeter continuum, where dust in the shell absorbs the 
FUV radiation emitted by the star cluster. 

To illustrate this, Figure 23 shows an IRAC band 4 
image highlighting the most massive of these clusters, 
LK11, with about 80 OB stars powering the DR22 region. 
Note that the PAH emission in this image is from the 
PDR, outside the ionization front, whereas the 24 /xm 
emission in Figure [12] shows greater correlation with the 
ionized gas. The BLAST sources occur along CO ridges 
to the north and south-east, the latter clearly arc-shaped 
(making this the "smiley nebula"). The most luminous 
source C92, coincident with a compact H II region, must 
be internally powered (see below). LK11 is shaping this 
clump and both the geometry and proximity suggest that 
star formation has been triggered there. 

Similar geometrical statements can be made about 
LK09 and C30 (and C27, C33) in DR7 (Fig. [1, 
LK12 and C115 in DR17 (Fig. |9k, LK13 and C127 m 



G81. 445+0.485 (the Diamond; Fig. 11), and DB22 and 
C182 (and C177, C178) in ECX6-33. 

Although these H II regions have evolved, the exciting 
clusters are still quite young; their stellar density is typi- 
cally five times higher than for evolved open clusters ( |Le| 
Duigou fc Knodlseder|[2002l ). Thus their radii R 50 (halt 



population radius), about V (0.5 pc), are perhaps indica- 
tive of the size of the precurs or molecular clump. The 
stellar masses are estimated by |Le Duigou fc Kn odlseder 
(|2002D to be 300 to 2000 M© (a few times more for 
JJR22) and so the precursor clump would be more mas- 
sive still. BLAST clumps are of this angular size, but 
not nearly as massive, not even up to these stellar masses 
let alone allowing for inefficiency of star formation. We 
conclude that there are presently no BLAST clumps ca- 
pable of forming such massive clusters. The embedding 
CO clumps are perhaps massive enough, but would have 
to condense considerably to form a concentrated cluster. 
The lifetime of these condensations could be quite short, 
lowering the likelihood of detecting this stage in a single 
GMC. 

5.3. Triggered Star Formation 

Spatially, a massive star cluster can promote further 
star formation in nearby molecular material by driv- 
ing winds and shocks, which both sweep up material 
and overrun and compress pre-existing condensations. 
Subsequent gravitational instability and/or radiatively- 
driven implosion of the overrun condensations collec- 
tively lead to what is known as "triggered" star formation 
( Elm egreen|1998| |Zinnecker fc Yorke||2007| |Koenig et aT] 
|2008p . Thus a sequence ol generations ol star formation 
can occur with spatial ordering (re cognized as "sequen- 
tial" star formation). For example, |Koenig et al. (2008) 
studied these phenomena in W5 with Spitzer. 

As discussed, in the vicinity of the OB star clusters in 
Cyg X there is evidence of further star formation. The 
most persuasive, if circumstantial, evidence for triggering 
is in the DR22, DR8, DR17, and DR15 molecular clouds. 

5.4. Compact H II Regions 

Earlier in their development, when massive protostars 
first become hot enough to emit ionizing UV photons, 
they initially ionize only a dense core, producing a hy- 
percompact H II region, potentially optically thick at 
lower radio frequencies. Subsequent expansion of the 
ionized gas leads successively to ultracompact, com pact, 
and then evolved H II regions ( |Churchwell 2002). The 
initial stage, which occurs when accretion is still strong, 
is accompanied by energetic bipolar outflows. 

Circumstellar dust reprocesses the absorbed radiation 
from the protostar, re-radiating the energy at infrared 
wavelengths. Thus in the Galactic Plane, embedded H II 
regions are the most luminous objects observed in the 
mid and far-infrared, first by IRAS and then MSX. 

Hot cores, dense mole cular clouds with te mperature 
> 100 K as classified by Kurtz et al.| ( |2000[ ), are often 
thought to be precursors ol UC H II regions. In our sur- 
vey, we do not detect any source warmer than ~ 40 K. 
Nevertheless, some of the brightest sources detected by 
BLAST in Cyg X are found among well-known H II re- 
gions, with previous studies revealing them to be com- 
pact or ultracompact; they are all detectable in 21-cm 
radio continuum emission in the CGPS. Five of these 



studied by |Motte et al. (2007) are tablulat ed in Table [5 ] 
and discussed in the next subsection. Motte et al. (200?) 
obtained masses of these sources assuming a dust tem- 
perature of 40 K. We are actually able to measure the 
temperature, finding 30 to 40 K. Still, it should be re- 
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Fig. 24. — Same as Fig. |17| but highlighting physical properties 
deduced from ancillary data. Many of the most luminous sources 
are ionizing and have associated radio continuum emission. Close 
to the nucl ear burning locus, these contain deeply embedded mas- 
sive stars (Motte et al. 2007): A HCH II; O UCH II; V CH II. 
Grey filled circles are other BLAST sources with compact 21-cm 
continuum emission, outside the survey area of MAMBO. Proto- 
stars in stage A are at somewhat smaller L/M than the c ompact 
H II r egions , near the accretion locus. These are classified by Motte 
et al. (2007) as □ HLIRPC; o MIRQP; and * IRQP. At the lowest 
L/M, clumps are externally heated (stage E), while at intermediate 
values the clumps are probably powered by low ma ss protostars; 
the thick dashed lines are the model predictions of Krumh olz fc| 
|McKee| (2008) (eq. [7] in § |OJ for surface densities = l.U, 0.1, 
and 0.01 g cm -2 . The vertical arrow signifies the result of evo- 
lution in L as star formation progresses in a clump. This might 
be interpr eted as the result o f the evolution of a single massive 
protostar (Moli nari et al.||2008| ). However, the initial rise s eem s 
more likely to be the result ofpower from low-mass YS Os (§ |6.2[). 
According to the model of |Krumholz McK ee (2008), unless Li 
exceeds a critical value near 1.0 g cm~ 2 , a ma ssive star will not 
form and so the rise in L will be arrested (§ |6.4[ ); this possibility is 
indicated schematically by the short horizontal arrow. Finally, the 
arrow segment pointing to the upper right represents qualit atively 



arrow segment pointing to the upper right represents qualitatively 
the evolution in the simulations by Smith et al. (2009) and Wang 
|et al.| (|2010), where mass is fed into gravitational potential minima 



trom larger scale structure. 

membered that BLAST measures dust emission includ- 
ing the outer envelope of these sources, which is cooler 
than dust near the core; with higher angular resolution 
MAMBO should sample slightly warmer dust on average. 

These five sources are marked in Figure 24 with special 
symbols. They are clearly internally powered, beyond 
the transition to predominant nuclear burning. For ref- 
erence, a 20 M zero-age mai n sequence star of spectra l 
type 09V has L bol = 10 5 L (|Schaerer & de Koter|l997|). 



mass of these clumps, > 100 M , they also seem likely 
to host more than a single star. 

The precursors to these sources should be cooler, at 
similar mass. In the range of a few 100 M , there are 
many such BLAST sources in Figure [24} the earlier stages 
to be discussed below. However, there do not seem to be 
the cooler equivalents of most massive sources (W75N, 
DR21, DR21 OH), suggesting that the precursor stage to 
the most massive (dense) clumps is relatively short-lived 
or that when cooler the clump is to o ext ended to have 
been classified as a BLAST source (§ |6.4). 



Also recorded in Tableware another 11 instances (in- 
cluding two distant ones) of BLAST sources with appar- 
ently ass ociated compact 21- cm radio emission, not tar- 
geted by |Motte et al.| ( |2007| ). These are slightly cooler, 
with median temperature 27 K. These are marked in Fig- 
ures [18] and 24 with grey filled circles. 

A histogram of the bolometric luminosities of all of 
these sources, including the luminosity corrections for 
the two distant ones, is given in Figure [25] The lumi- 
nosities of the five marked sources are atthe high end, 
whereas most of the others have much lower luminosities. 
The luminosities can be used to calculate the correspond- 
ing photoionizing flux Qo and check for c onsi stent radio 
emission (assuming no optical depth; see §|3.2|). The low- 



est luminosity sources in Figure 25 would not account for 
the associated H II emission. Inspection of the IRAC im- 
age reveals all but C76 to have comet ary structure (those 
of higher mass and luminosity were already mentioned 
above: C55, C73, C80; all but C14 appear to be influ- 
enced by Cyg OB2) and our interpretation is that the 
radio emission comes from external ionization. C76 lies 
in close projection to the 07 V star MT771. 

5.4.1. Specific Embedded H II Regions in Cyg X 

Here we make brief commen ts on the five marked 
sources from Motte et al. (2007) in Table [5j plus three 
others that appear to be m the same advanced stage of 
evolution. 

CI, AFGL 2591: CI coincides with AFGL 2591, an 
UC H II r egion that is among the most luminous sources 
in Cyg X. |Motte et all] fl2007| identify a core S26 within 
a clump So! We rind dust temperature 40 K, b olomet 



ric lumin osity 2.7 x 10 L , and mass 190 M . Camp 



bell| Ql984|) discovered that AFGL 2591 is comprised oi a 



young stellar group and deduced that the H II region is 
being generate d by B0 stars. This region also has pow - 

^19921). 



erful outflows (Bally & Lada 1983; [Poetzel et aL 
Schneider et al. Q2006 ) argue that the cometary s 



Even the most luminous of these embedded stars at the 
distance of Cyg OB2 is not that luminous. Given the 



nape is 

due to the influence of Cyg OB 2. 

C88, DR20: Compared to CI, C88 is slightly cooler 
(32 K), less luminous (8.6 x 10 3 L ), but equally mas- 
sive (190 Mf? ) ). It c orresponds to clump N4 (core N10) of 
|Motte et al. (2007). The embedding region (see Fig. 10) 
is bright in 21-cm radio continuum, mid-infrared, and 
PAH emission. C88 and C87 are two bright sources along 
a ridge containing C90. C88 is closest to the CGPS radio 
peak, but offset to lower longit ude, whereas the U C H II 
region G80.86+0.4 observed by |Kurtz et al.| ( p^994| ) is cen- 
tered on C88. C8 8 also c oincides with 1KAS 20350+4126. 
|Odenwald et al. (1990) suggest excitation by an 05-6 
ZAMS star, but this seems too early. 

The source C87 (clump N6, core N14) is relatively cold, 
27 K, with no free-free emission. Its position on the L- 
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Fig. 25. — Histogram of Lboi ^ or BLAST sources with compact 
H II emission. Two of the most luminous (st ripes) are at a larger 
distance than Cyg OB2. The five noted by |Motte et all (2007) 
are shaded black. Ot her lower luminosity sources appear to have 
external ionization (§ \5A\. 



M plot suggests most of the luminosity is generated from 
accretion (see below). 

C92, DR22: C92, also known as IRAS 20375+4109, is 
a co mpact H II region in DR22 with a massive dense core, 
N58 flMotte et al.|2007| ). We find dust temperature 39 K, 
bolometric luminosity 2.4 x 10 4 L , and mass 220 M . 
Odenwald et al. (1986|) studied this region at both radio 



and tar- infrared wavelengths and assuming it to be 45 K 
and 3 kpc distant, they deduced that a single 06 ZAMS 
star was accounting for a luminosity of 3 x 10 5 L . This 
would be relaxed for the closer distance adopted here, 
and as mentioned above th ere i s likely more than a single 
star in a BLAST clump (§[5^4). 

CI 55, DR21: CI 55 is the most luminous source in the 
BLAST survey region with dust temperature 36 K, lu- 
minosity 7.0 x 10 4 L , and mass 880 M . DR21 has a 
dense core N46 with signatures of outflows (iMotte et al.l 



2007), suggestive o f an early stage OB cluster. Support 



lve of multiplicity, |Motte et al.| ( 2007 ) detected another 
two cores within a r radius in tneir clump N15; our mass 
agrees with the total mass of the three cores. This clump 
in the DR21 cloud complex (see Fig. 11) harbors one of 



the most studie d UC H II regions: microw ave ( jDownes 
|Rinehart"l|1966|); ra dio ( |Kurtz et al.||1994[); near- i nfrared 



Schneider et al. 



(2006). 



payis et al.||2007[ ); and L6 C0 ^ 

Stud ying hydrogen recombin ation and the ammonium 
line, [Cyganowski et al" ( 2003 ) found cometary morphol- 
ogy associated with the H II region. It is also identified 



as MSX6C G81.6802+0.5405, an MSX source, specifi- 
cally bright at 8 /im, and has saturated Spitzer so urces. 
There is an embedded near-infrared cluster, W75S (iBical 



et al. 2003). The position of C155 in the L-M diagram 
confirms that star formation is well established. 

CI 69, W75N: The W75N cloud complex shown in Fig- 
ure [TT] is a well-known massive star forming region. Af- 
ter "Ers^_being detect e d at low resolution in the radio 
by |Westerhout| fll958|); |Downes fc Rinehart| (p66b, this 



complex was studied to decipher numerous protostellar 
physical conditions." A wide variety of sources was found 
embedded in this c l oud, e.g., a group of infrared source s 

2006} [ Davis ■ et al. | |2007[) 



(Moore et al. 1991 Persi et al. 



UC H II regions, and H^O an d OLL niasers ( [Hunter et al 



1994). |Haschick et al. (1981) identified three ionized re- 
gI5ns within W751N, namely W75N (A); W75N (B); and 



W75N (C). Subsequent higher resolution |Hunter et al 
([1994[) resolved W7 5N (B) into another three subregions 
Ivlotte et al.| ( |2007 ) found three sources associated with 
the W75N cJoTlGT^lump N13), of which W75N (B) is the 
most massive object, with a dense core N30. They clas- 
sified this source as a HC H II region. BLAST with V 
resolution observes only one source. The derived proper- 
ties are dust temperature 36 K, luminosity 6.2 x 10 4 L , 
and mass 760 M , very close to C155. The radio emis- 
sion at 21-cm is however much less, probably because of 
self absorption (there i s appreciable self-abso rption in the 
spectrum of DR21 too ,|Wendker et al.|1991| ) but possibly 
because of a different IMF in the embedded clusters. 

There is another embedded cluster DB20 a few arc 
minutes to the north, coincident with clump N14 of 



Motte et aU ( )2007P and its core N28. 

U32, DKlb cloud: The DR15 co mplex w a s prev iously 
observed at rad i o wa velengths by Colley (1980), and 
Odenwald et al. ( 1990[ ), who predicted that it hosts OB 
stars. Both tree- tree emission and PAH emission ob- 
served by MSX suggest on-going star formation. The 
morphology of this cloud is complex, with extended 
emission on a range of scales. We could resolve two 
sources, C32 and C28. The cluster LK08 a nd ra- 
dio peak ar e located between these (Fig. [7]). |Motte| 
et al.| (|2007 ) detected three embedded sources, of which 
two7MSXT9.2963+0.2835 and MSX 79.3070+0.2768 are 
within C32. The Spitzer sources are saturated. We find 
dust temperat ure 46 K, luminosit y 3.3 x 10 4 L , and 
mass 110 M . |Motte et al] ( |2QQ7j ) designated core S41 
in clump Sll (corresponding to C32) as a HLIRPC. Our 
analysis based on the location of C32 in the L-M plot 
suggests a slightly later stage like the other embedded 
H II regions being discussed. However, it is a complex 
region that will benefit from better resolution as provided 
by Herschel. It is possible that the temperature of the 
dust is influenced by the nearby cluster LK08, and fitting 
a single temperature SED for this source would probably 
underestimate the clump mass. 

CIS, DR 18: DR18 was studied exte nsively in the near 
infrared by Comeron & Torra (1999) who found a stel- 
lar group. Most of the illumination is due to a single 
B0. 5 V star which can be observed in the visible as well 
as the infrared. The CGPS radio peak is at the head 
of the larger cometary structure (see above) and proba- 
bly is affected by external ionization. Consistent with the 
cometary structure, C73 lies between the ionized gas and 
the CO clump. C73 has properties very similar to C88 
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in DR20. This intere sting region is outside the coverage 
of |Motte et aT| J2QQ7b- 

C3U, JJK7: Figure 6] shows LK09 located in the centre 
of the DR7 H II region arc. The CGPS H II emission 
peak and ECX6-18 lie between BLAST sources C30 and 
C27, the f ormer much more m assive. This region is not 
covered by Motte et al. (2007). The derived parameters 
for C30 are very close to those for C32 , as suming the 
same distance. However, as remarked in § |3.6[ its velocity 
places it well behind the Cyg O B2 complex, at 3.6 to 
7.5 kpc (we adopted 6.8 kpc in 



4.2). Thus this source 



should be about 16 times more massive and luminous, 
as listed in Table [5j ranking it with the most powerful 
sources in this survey. 

5.5. Stage A: Accretion Power, Protostars, and MYSOs 

As mentioned, BLAST cannot hope to see single stars, 
though the clump might be dominated by the most mas- 
sive member of the group. During the next earlier stage, 
precursor to the embedded H II regions, matter collapses 
under gravity at multiple sit es inside the clump, which 
are newly born stellar cores ( |Mac Low fc K lessen 2004). 
The luminosity of the clump originates from the cumula- 
tive energy of accretion onto cores. During this accretion 
stage, protostars build up luminosity with little loss of 
envelope mass. In our L-M diagram (Fig. 24), the evolu- 



tion should be vertical, crossing the regionoetween the 
dotted red line (L/M = 1 L /M ) and the dashed green 
line. The former is coincident ally close to diffuse cirrus 
emissio n in equilibrium wi th interstellar radiation field at 
17.5 K flBoulanger et al.||1996|). The latter, taken from 



Figure 9 of |Molinari et aL| ( |2008 ), is a best fit to class 



sources. Sources plotting in this region would have cooler 
dust temperatures, in the range 18 — 25 K and L/M 
> 1 Lq/Mq would clearly indicate internally powered 
sources. 



Motte et al. (2007) have identified several massive in- 
frared quiet protostellar cores (MIRQP), sources with 
mass greater than 40 M with at least some evidence 
of stellar activity and yet weak emission in the mid in- 
frared. All of the massive infrared-quiet cores are associ- 
ated with high velocity SiO emission, providing convinc- 
ing evidence of powerful outflows and on-going accretion. 
The SiO emission of MIRQPs is typically stronger than 
from high luminosity infrared protostellar cores (HLIR- 
PCs), indicating that the former are in a more vigorous 
accretion phase. 

Of the 17 MIRQP, 8 are within our BLAST survey area 
and are readily detected and characterized (Table pj. In 
the L-M evolutionary diagram, they occupy the pre- 
dicted region (see diamonds in Fig. 24). Judging from 
Figure |24j there are many other BLAdT sources with 
similar properties. The sensitive Spitzer images reveal 
the sub- structure, with many accreting cores over a range 
of fluxes (see also 



6.2). 



CI 57, DR 21 (OHJT The most luminous MIRQP is 
in DR21(OH), which is lacking in 21-cm radio free- free 
emission. However, high-sensitivity subarcsecond ob- 
servations w i th th e VLA at centimeter wavelengths by 
Araya et al. (2009) reveal a cluster of radio sources; the 



strongest emission is toward the molecular core MM1, 
but is perhaps from shock- ionized gas in a jet rather than 
from a nascent compact H II region. C157 is also cooler 
than the strongly ionizing sources discussed above, as 



one might expect. It has a dense core N14 a long with 
five neighboring fragment s ([Motte et al.||2007 ) in clump 
N14. |Harvey-Smith et al.| ( |zuU8| ) have observed 6.7 GHz 
methanol emission in all five spots, but only two of them 
(DR21(OH) and DR21(OH)N) exhibited strong peaks. 
Our limited resolution integrate s the dust em i ssion from 
all nearby fragments. Though |Motte et al. (2007) has 
classified the core as MIRQP, we note that the signif- 
icant bolometric luminosity of C157, 2.9 x 10 4 L , is 
probably indicative of the entire cluster, detected in the 
near-infrared as DB19. 

The next most luminous MIRQPs are in C7 and CI 15. 
C7 is an isolated source whereas CI 15 is at the tip of 
the DR17 molecular pillar. CI 19 is similarly at the tip 
of an adjacent pillar (Fig. [9]). C36, among the coolest, 
is just north of DR15, at the most e volved end of the 
IRDC G79.34+0.33 discussed below(§ [5A 2|. There is a 
variet y of indications of star formation ( jRedman et aL 
2003[), and Spitzer images indicate a rich stellar group. 
There are two sources classified as infrare d quiet pro- 



tostar (IRQP, not so massive as MIRQP) by Motte et al. 
(2007) which appear to BLAST to be as luminous as the 
MIRQP- containing sources. The more luminous is C104 
(N62 in clump N20 to the north-east of DR22) in the 
molecular filament (velocity —5 kms -1 ) extending from 
DR22 through MIRQP C107 to MIRQPs C141, C145 and 
cluster LK15 in DR23 (see Fig. 12). The most prominent 
molecular clump along this filament contains CI 38 (N69, 
in clump N23); it defines the boundary of brightest ion- 
ized part of DR23, powered by cluster DB17 to the east 
(Fig. [12). The other is C167 (N24 in clump N10 to the 
north-west of W75N). At CO velocity —3 kms -1 , this is 
in the DR21 cloud complex (Fig. [Tl|). The Spitzer im- 
ages indicate widespread embeddeclstar formation in the 
C167 cloud and adjacent BLAST sources. 

The two HLIRPC (other than C32 discussed above) 
plot toward the high luminosity (high temperature) side 
of this region in the L-M plane. These correspond to 
the BLAST sources C85 and C87. C85 contains source 
N6 in clump N3 of [ Motte et all (|2007D . It coincides 
with DR20NW Schneider et al.||2UU6| Tat the head of 
the above-mentioned elephant trunk. There is an ion- 
ized rim on the face toward Cyg OB2, but apparently 
little internal ionization. Its position on the L-M plot 
suggests most of the luminosity is generated from ac- 
cretion. Compared to C88 also in the DR20 molecular 
cloud, C87 (clump N6, core N14) is relatively cold, 27 K, 
with no free- free emission. 

5.6. Stage E: Externally-heated Cold Early Stage 

Prior to the stage of active accretion, there must pre- 
exist a cold molecular clump with only external heat- 
ing. A low temperature and corresponding low lumi- 
nosity to mass ratio (L/M < 1 L /M ) would be 
diagnostic of this stage. As a forerunner of SPIRE on 
Herschel, BLAST was designed to search for the submil- 
limeter emission from such objects in the earliest stages 
of pre-stellar evolution, and according to Figure [17] has 
been quite successful, with 28 sources down to the ap- 
proximate 250 jam completeness line of about 30 Jy. 
These compact sources are sufficiently cold and massive 
to be above the nominal Bonnor-Ebert critical mass. The 
higher mass objects can be considered as potential pre- 
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Fig. 26. — Left: cool BLAST sources at 250 /im. Right: corre- 
sponding 24 am MIPS images, showing evidence for nascent star 
formation. Top: C116. Bottom: C36 at the transition to MIRQP, 
in IRDC G79.34+0.33. Radius of circles 1.3'. 

cursors of the MIRQP accretion-powered stage. 

5.6.1. Cold Clouds with Nascent Star Formation 

It is a general consensus that massive stars should be- 
gin to form within stellar nurseries deeply embedded in 
a dense envelope of dust and molecular gas. Because 
Spitzer is so sensitive and probes at wavelengths where 
the dust can be penetrated, evidence for such nascent 
star formation (not a lot of total power) is readily found. 
Figure 26 shows two examples. One is a cold BLAST 
source C116 (alon g the DR22 to DR2 3 filament) cor- 
responding to the |Motte et al | ( )2007| ) core N70. We 
find dust temperature 17.1±0.8 Iv, bolometric luminos- 
ity (2.2 ± 0.4) x 10 2 L , and mass (1.6 ± 0.4) x 10 2 M . 
The other is C 36, already containing a lo w tem perature 
MIRQP (§ [575]), embedded in an IRDC (§ [5A2| . 

5.6.2. IRDCs 

Given the appropriate geometry, highly dense, massive 
molecular clouds with significantly large extinction can 
be seen in silhouette as infrared dark clouds (IRDC). 
With the spatial resolution and sensitivity of MSX, many 
IRDC s have been found at Band A (8.6 /im) ( |Simon et al.| 
|2006[ ). The signature of an IRDC without a developed 
protostar is l ack of dust emiss ion at mid-infrared IRAS 
wavelengths (Egan et al. 1998|. However, that does not 



mean there is no low- power nascent star formation, as 
illustrated below. 

Because of the high column density, BLAST sees 
prominent IRDCs in emission. An example is the ridge 

which contains the IRDCs 



(filament) shown in Figure 27 

G79.34+0.33 and G79.27+(T38 ( |Egan et al.|[l998| > with 



corresponding BLAST clumps C36 and C26, respectively. 
Although it was initially suspected that G 79.34+0.33 was 
assoc iated with the DR15 cloud, later |Redman et al.| 
(2003) found no such evidence. This filamentary ridge 
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Fig. 27. — Multi-wavelength views of the IRDC ridge in Cyg X 
near I = 79°. 23 and b = 0°.45. Each thumbnail has dimension 
12'x 16'. 

is associated with distinct 13 CO emission at velocities 
near 1 kms" 1 (Fig. [7]). G7 9.27+0.38 in the ce ntral part 
of the ridge was studied by Wu fc Yang| ( |2005[ ) . 

Finer details can be found with the Spitzer images, 
both in extinction and in the presence of point sources. 
Both of these BLAST sources show evidence of nascent 
i and ! 



star formation (Figs. 26 and 27). In fact there are embed- 
ded mid-infrared objects right along the filament, even 
the upper part of the ridge where the less massive con- 
densations are C25, C24, and C23. 

Note that this ridge is seen in e xtinction against th e 
shell of the LBV star G79.29+0.46 ( |Wendker et al.|1991[ ), 
providing relative distance information which locates the 
star behind the Cyg X region. The shell, even the faint 
outer halo seen at 24 /xm, shows no signs of interaction 
with Cyg OB 2, ind icating considerable separation. |Krae- 



mer et al. (2010b) have quantified this using the stellar 
parameters, obtaining a distance of 3 kpc. 

Figure [28] shows an interesting long filament that shows 
up well in the 24 /im image, in emission but changing to 
absorption in front of a bright band of cirrus. At IRAC 
8 /im it is only in emission, the bright cirrus band being 
absent. This structure is detected in BLAST emission 
as well as CO, extending the DR21 ridge to the south 
(Fig. 11). The characteristic CO velocity, —3 kms -1 , 



confirms that this is related to the DR21 molecular cloud. 
The extended mid-infrared source in the upper part is 
apparently warm and has no BLAST counterpart. 

5.6.3. Starless Cold Clouds 
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Fig. 28. — MIPS 24 /im image showing dark lane exte ndin g from 
DR21 ridge. BLAST sees t his lane in emi ssion (Fig. [TT). The 
contours are 13 CO emission ( jSchneider et al.|2006) integrated over 
-7 to 1 kms" 1 . 

Focusing on the earliest stage, are there any sources 
that are so far starless? The best candidates to search 
should be those in the region with L/M < 1 L /M . 
In terms of internally generated power, they all can be 
considered "starless." But when examined in the MIPS 
24 /am image, most of these candidates do in fact have 
point sources, indicating nascent sta r formation though 
probably only low mass YSOs (§ |6.2[ ). 

However, we found one instance that is apparently 
starless at 24 /im, C81, an isolated feeble source with a 
weak CO signature at +6 kms -1 . While there are some 
point sources at IRAC 8 /im, these do not have spectra 
rising to 24 /mi as might be expected of YSOs; also there 
is no special concentration to the BLAST source and so 
they are probably field star contamination. The IRAC 
image also shows several small IRDCs. 

Two other possibilities are CI 03 and C121. (Curiously, 
C103 exhibits diffuse emission at 24 jam including being 
crossed by a partial thin shell of radi us 2.8' centered at 
G81.0425-0.1155.) C108, clump N2 of |Motte et aL| p007|j 
with no core, is crossed by IRDCs at 24 and 8 JIm~. and 
might be another good example of a starless clump. 

6. DISCUSSION 



The L-M diagram (Figs. 17 and 24) can be adopted 
as a diagnostic tool for describing the early evolutionary 
stages of massive star formation (Molinari et al.|[2QQ8[). 



With our sample of compact sources now covering a 
broad range of luminosity and mass, particularly extend- 
ing to low L/M ', we are in a position to further assess its 
utility. We first review how independent empirical evi- 



dence for evolutionary stages relates to the characteristic 
position in the L-M diagram and then, to the extent pos- 
sible, use theoretical models to provide insight into the 
evolution. 

6.1. High L/M 

In the Cyg X region, 16 BLAST sources have de- 
tectable compact 21-cm continuum emission. All have 
strong 100 /mi emission and so are not background radio 
galaxies. Our detections of these sources are represented 
by filled light-grey circles in Figure [24] Five of these 
were al so noted and class ified for H II region compact- 
ness by |Motte et al.| ( |2007| ) and are marked as triangles, 
circles, and inverted triangles. All of these clumps ap- 
pear in the higher T (> 30 K) and L/M region, near 
the nuclear burning locus, above the stage A sources and 
the empirically-obtained accretion line. An extensive dis- 
cussion of individual e mbe dded H II regions in Cyg X 
has been presented in § 5.4 High mass protostars inside 



these massive clumps account for the high ionization rate 
and luminosity. We do not have the angular resolution 
to distinguish multiplicity, but the clumps contain multi- 
pl e sou rces when examined in the near and mid infrared 

(r 



6.2). 
5KJ (2 



1) emission is often used as a tracer of 
shocked gas, associated with outflows which are an in- 
trinsic consequence of accretion process. After the onset 
of nuclear burning stage, massive stars build up an in- 
tense radiation pressure which halts further accretion. In 
accordance with this expectation, the luminous sources 
DR21 an d W75N have only low intensity SiO (2 1) 
emission ( [Motte et al.||2007[ ). In fact, all of their high 
luminosity 1R Protostellar cores (HLIRPC) correspond 
to BLAST sources with high L/M and also do not have 
a striking molecular emission accretion signature. 

Slightly earlier in the evolution, clumps are already 
quite luminous, near the empirical accretion locus for 
massive stars; single-star cores within them would be 
classified as class 0, at least for lower mass stars. On 
the basis of their position in this region in the L-M plot, 
we would expect these BLAST sources to show indepen- 
dent and direct evi denc e of (bipolar) outflow and indeed 
this is the case (§ |5.5| ). |Motte et al.| ( |2007| ) identified 
a class of massive lR-quiet protostellar cores (MIRQP) 
with strong SiO (2— >4) emission, indicative of on-going 
active accretion. The eight BLAST sources correspond- 
ing to their MIRQP are marked by diamonds in Figure 24 
and all of them lie in the stage A accretion-dominatecT 
region of the diagram. 

6.2. Low L/M 

In the Cyg X region we have discovered a significant 
fraction of the clumps with low L/M devoid of massive 
protostellar cores. This low L/M regime is for us the 
most intriguing part of the L-M diagram, potentially 
containing th e story of the ear liest stages of massive star 
formation. |Molinari et al.| ( |2008[ ) show vertical tracks cor- 
responding to~tne^voIutionoTa single and massive pro- 
tost ar, and this is suggested by the vertical arrow in our 
L-M diagram. However, this is probably a misleading 
interpretation for low L/M. Even though these clumps 
appear above the L L = M' line, and are arguably self- 
luminous, there is no obvious signature of a high mass 
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Fig. 29. Number of YSOs (classes 0, I, and II) per 100 M© 
vs. L/M. The cross, diamond, and plus symbols represent surface 
densities E in the ranges 0.005 to 0.02, 0.02 to 0.08, and 0.08 to 
0.4 g cm -2 , respectively. The vertical (dashed) lin e corresponds 
to 20 K. a ty pical value where E in the model o f |Krumholz Sz\ 
McKe e ( 2 008) is about 1 g cm -2 for the masses of BLAST sources 
(see §|Qfancl Fig. [24} . 

star being formed. Given that massive protostars prob- 
ably evolve rapidly through this region, it is statistically 
improbable to observe this stage. 

The more likely alternative energy source is the col- 
lective power of many lower mass YSOs. In § |5.6.1 



have already previewed some evidence for multiple YSOs 
in cold massive clumps. Since the BLAST observations 
were carried out, a c ensus of YSOs has been obtained by 
Beerer et al. ( 2Q1Q[ ) based on SEDs from Spitzer imag- 
ing of the northern part of the Cyg X region. Using 
their map of Spitzer YSO positions (their Fig. 8) we find 
YSOs, often multiple YSOs, associated with most of our 
clumps. For the brighter sources, we probably underes- 
timate the number of (crowded) YSOs when using their 
Figure 8 rather than a catalog. Figure [29] demonstrates 
that the observed L/M is correlated with the total num- 
ber of YSOs per unit clump mass, as expected energeti- 
cally. This also suggests considerable fragmentation into 
low-mass stars, an inevitable result of the ini tially low 
Jeans' masses in the cold gas ( |KrumhoTz||2QQ6[ ) . 

There are (at least) two hypotheses tor massive star 
formation, namely the competitive accretion model 
( Bonnell et al.|1998| ) where massive stars are formed due 
to dynamical coalescence of low-mass stars in a dense 
clustered med ium, and the turbulent core model ( McKee 
& Tan 2003) in which massive cores with non-thermal 



200 400 600 

BLAST Clump Mass [M ] 

Fig. 30. — Sum of the masses of cores within the BLAST clumps 
versus the clump mass. The linear correlation shown has a slope 
0.20. 

static way. Both of these theories predict that high mass 
star formation is preceded by an epoch of the low-mass 
star formation, converting a larger proportion of the to- 
tal mass in the competitive accretion model. Owing to 
low resolution, BLAST observations cannot differentiate 
between these two hypotheses on this basis. Nor do we 
have masses for the Spitzer YSOs. However, the inter- 
mediate resolution MAMBO data resolves sub-structures 
(cores) of size ~ 0.1 pc within our clumps. In Figure 30 
we compare the total mass of the cores within each clump 
to the clump mass (assuming the cores are at the same 
temperature as BLAST observes for the clump). We find 
that (only) 20% of the clump mass is in cores, however, 
four bright sources (W75N, DR21, DR21(OH) and C160) 
were excluded in the fit. 

6.3. Relationship to Surface Density 



According to the model of Krumholz & McKee (2008), 
for the low L/M regime the steady-state accretion lumi- 
nosity generated by the low-mass YSOs is: 



390 



M 



1 g cm- 2 100 M 



0.67 



L^ 



(7) 



support provided by the turbulence evolve in a quasi- 



In order to show the relative position of BLAST sources 
with respect to this model prediction, we plot loci in the 
L-M diagram of Figure [24] for £ = 1.0, 0.1, and 0.01 
g cm -2 . 

If this theory is adopted then we can directly deter- 
mine surface density from equation ^ for any source in 
the L-M diagram. This is the surface density that would 
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Fig. 31. — Observed surface density of BLAST clumps with 
L/M < 6 Lq/Mq vs. the su rface densit y requ ired in the low- 
mass YSOs model of |Krumholz fc McKee| (|2008 ) to produce the 
observed L given the observed M . Sources above the one-to-one 
line (dashed) are under-luminous compared to the prediction from 
the model, and so possibly very young. The source shown by a 
cross is near the map edge at 500 |Um. Sources with model surface 
density more than 1 g cm -2 are highlighted. 

be required for a clump of mass M to produce a lumi- 
nosity of L. However, we can also measure the surface 
density directly. In Figure 3lJ we plot surface densities 
of BLAST sources versus the model surface density, for 
sources below L/M < 6 L /M where most of the lu- 
minosity arguably comes from low-mass stars. The mea- 
sured surface density of the most of the sources is below 
the one-to-one line. This might be because we have un- 
derestimated E owing to poor angular resolution; qual- 
itative support for this is the evidence of sub-structure 
(§ 6.2). However, there are about 8 sources with surface 



density above the one-to-one line. These are very inter- 
esting, possibly a very early stage where the presumed 
steady state of low mass YSO formation has not been 
achieved. We have confirmed that there is a lower num- 
ber of YSOs associated with these sources, except for the 
source shown by a cross in Figure 31 which was difficult 
t o measure near the edge of th e image at 500 /im. 

Krumholz & McKee (2008) note that because cores 
with higher E are warmer due to the power of the low- 
mass YSOs, there is a critical E such that by the time 
the steady-state heating is achieved fragmentation into 
further low-mass YSOs halts. The threshold surface den- 
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sity is about 1 g cm -2 and this is plotted in Figure 
This characteristic surface density should be a marker 
of (potential) massive star forming regions, and there is 



some evidence consistent with this, at least for evolved 
objects. In the BLAST survey we find that the few most 
luminous sources containing compact H II regions have 
such a high surface density. The MAMBO cores contain- 
ing compact H II regions and MIRQPs with strong SiO 
emission also hav e a high measured su rface density (Ta- 
ble [3]). Similarly, Garay et al. (2006) found an average 



surface density of 0.8 g cm -2 for dense cores of ultra- 
compact H II regions. We find that the typical surface 
density of massive embedded st ar clusters observed by 
|Le Duigou fc Knodlseder| ( |2002[ ) in the Cyg X region is 
about 1 g cm -2 (allowing for the presence of some addi- 
tional gas - now dispersed - in the precursor). 

An observed high E associated with MYSOs might not 
be a signature unique to the turbulent co re model, or 
even a prerequisite. In the simulations by |Smith et al. 
(2009) of a competitive accretion model, "massive stars 



were not formed from a single massive thermally sup- 
ported fragment, but instead from a smaller core which 
accreted additional material channelled towards it by the 
potential of th e forming ste l lar clu ster." Likewise, in the 
simulations by |Wang et al.| (l2010| ) of a model of "clump- 
fed massive star formation regulated by outflow feed- 
back (ORCF), "the most massive object is not formed 
out of a pre-existing dense core" but rathe r "is controlled 
by the global clump dynamics." Recently Csen geri et al. 
2010[) observed a few of the massive young dense cores 
TUT, C115, C119, DR21(OH), DR21) in Cyg X to 
study the kinematic properties of the underlying envi- 
ronment. They found that non-thermal support due to 
turbulence at the scale of the dense core is not enough 
for a quasi-static evolution of the massiv e protostars as 
predic ted by the turbulent core model of |McKee fc Tan| 
(2003), and they favor competitive accretion. 

Therefore, a key question in investigating the mode of 
evolution is whether there is evidence for cold high-E 
clumps or cores before high mass protostar(s), and even 
the preceding low-mass YSOs, are established. Even 
if the evolution timescale were favorably long, this is 
difficult with BLAST, as can be seen directly from 
equation (|3|: at the spatial resolution of 0.5 pc about 
4000 Mq is needed to produce a surface density of 
1 g cm . We have noted that we do not see such mas- 
sive cold objects; if such dense cold objects do exist at 
lower mass, we cannot identify them by directly measur- 
ing their true E. Higher resolution is required for such de- 
tections; for example, MAMBO observations at the reso- 
lution of 0.1 pc in Cyg X need a mass of only 150 M to 
produce equivalent contra st. The predomin antly low sur- 
face densities obtained by |Elia et ah] ( |2010[ ) suggest that 
even the Herschel resolution is not sufficient for finding 
such dense sub-structures in most cold sources. 



6.4. Vertical Evolution in the L-M Diagram? 

Given the different theories, it is still unclear to us 
which of the low L/M BLAST clumps will give rise to 
massive star formation. If there is a universal IMF, then 
statistically a massive clump will be required to give birth 
to massive protostars. Therefore, the massive clumps 
identified by BLAST are certainly of interest in this re- 
gard. Naively, the evolution in the L-M diagram would 
be vertical (constant M) as the luminosity increases from 
accretion onto the massive protostar and then nuclear 
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burning. Even this is complicated because the power 
produced by low mass YSOs in most theories can mask 
the earliest low L stage of the massive protostar. 

However, it does not follow that massive clumps/cores 
always form massiv e star s. In fact, the model by 
Krumholz & McKee] (|2QQ8|) explicitly predicts that this 
is above a critical value of about 



will not happen unless . 
1 g cm -2 . Otherwise, their prediction is that the clump 
is destined to continue to produce only low mass stars, 
at some presumed steady-state rate. Thus the vertical 
rise in L is arrested at a ceiling for a given E. We have 
marked such a possibility in Figure 24] with a horizon- 
tal red arrow at a level which would correspond to the 
sub-critical E of the clump. The end result is a cluster 
or association of only low mass stars, i.e., one with a 
truncated IMF. Because the precursor is of low E the re- 
sulting assemblage of stars must be as well, making them 
more difficult to detect, and there is no high mass star to 
draw attention to the region by its high luminosity and 
ionization. 

The empirical context, a snapshot of the stages of evo- 
lution, is that there is a hierarchy of structures. BLAST 
identifies compact sources of size about .5 pc within the 
CO c lumps of spatial extent about 10 pc ( jSchneider et al.| 



2006). Within most of the BLAST clumps, MAMBO is 
able to find smaller and denser sub-structures of scale 
0.1 pc. Presumably these structures have been created 
and evolve by accumulating mass in the dynamical envi- 
ronment of the larger reservoir and possibly also simul- 
taneously transfer mass to the immediately lower scale 
sub- structures. At any given time these structures will 
each have a characteristic E. At some point does this 
become frozen in time so that we can discuss the fate of 
objects of fixed E? Or tracking a clump, could E and M 
evolve, allowing evolution toward the upper right in the 
L-M diagram? We have indicated this latter possibility 
of course only schematically. Re- 



as well in Figure 24 



cent n umerical simulation s of the competitive accretion 
model ([Smith et al.||2009[ ), and the ORCF model flWang| 
|et al.||2010| ) which has features of both competitive ac- 
cretion and the turbulent core model, show a channeling 
of mass to bound potential energy minima of the parsec- 
scale dense clumps studied, emphasizing the importance 
of clump- level dynamics. In these models, when the ma- 
terial is cold it is more distributed and so would not be 
identified as a high mass, high E, cold clump. By the 
time that mass builds up in massive cores, accreting on 
massive protostars, the compact source would have low 
mass YSOs as well and have been warmed up. Thus the 
cold precursors of dense clusters might be more extended 
objects than clumps, a massive reservoir yet to be chan- 
neled by gravit y into a higher surface density state. Some 
cold IRDCs (§ ]5lT2l Fig.|3|, within which BLAST iden- 
tifies substructure as sources, might fit into this scenario, 
but dynamical diagnostics would be needed to pursue the 
relevance of this possibility. 

7. CONCLUSION AND FUTURE WORK 

An unbiased survey conducted by BLAST has enabled 
us to detect both compact sources and the embedding 
diffuse emission as a basis for studying the earliest stages 
of star formation. We have quantified contrasting dust 
temperatures of diffuse extended structures by correlat- 
ing the fluctuating emission in IRAS and BLAST bands. 



To reveal the global morphology, we have studied the re- 
lationship between dust continuum emission and radio, 
mid-infrared, and 13 CO line emission. Comparing the 
total mass present in the BLAS T compact sources with 



the parent CO clump mass (Schneider et al. 2006) 



we have obtained a linear correlation with slope of 0.05, 
which might be indicative of the star formation efficiency 
but also relates to the relative lifetime of the compact 
clump phase. With mult i- wavelength photometry we 
have tightly constrained the dust emission SEDs of the 
compact sources to measure temperature, luminosity and 
mass. We have diagnosed the earliest stages of star for- 
mation in Cyg X, and the evolution, in the context of 
the luminosity - mass diagram of the compact sources. 
We have studied eight DR regions which are known as 
massive star formation sites in Cyg X. We have identified 
16 sources which have compact H II emission and found 
an average dust temperature of 33 K. Not all of these 
identified sources have sufficient luminosity to be ioniz- 
ing sources, pointing to the influence of Cyg OB2 creating 
external ionization. The complex morphology along with 
the direction of cometary shapes of some compact sources 
provide direct evidence of interaction with the massive 
Cyg OB2 association. Along the periphery of the evac- 
uated H II region created by some of the more evolved 
stellar clusters, there is a signature of triggered star for- 
mation. We have observed a complex of infrared quiet 
dark clouds where star formation has already started to 
take place. To BLAST, and also SCUBA and MAMBO, 
this appears as an impressive ridge of cold dust emission. 
However, in sensitive high resolution 24 /am MIPS and 
8 jam IRAC images, there is evidence of a stellar nurs- 
ery, indicating a potential site of massive star formation. 
In our search for the precursors of clusters, however, we 
have not found any cold compact clump which is massive 
enough to become a star cluster like the Cyg OB2 asso- 
ciation or even the other lesser embedded clusters. In 
the future, the unbiased Planck Cold Core survey might 
find such objects in unevolved GMCs, possibly in a more 
extended lower surface density state than in the final 
cluster. Resolving multi-scale structures is a new chal- 
lenge underlying investigations of the earliest stages of 
star formation. BLAST with its limited resolution could 
only resolve spatial structures of 0.5 pc at 1.7 kpc away. 
On the other hand, MAMBO with its 11" resolution was 
successful in separating sub- structures on a scale of 0.1 
pc; further progress should be possible with Herschel. 
To determine the role of surface density and more gener- 
ally which cold high mass clumps will evolve to produce 
high mass protostars even higher resolution observations 
will be required, and so there appears to be a rich fu- 
ture in uncovering th e unseen story using SCUBA2 (8" 
resolution a t 450 //m, |Holland et atl |2006), CCAT (3.5" 
at 350 /mi, |Sebring||2010[ ), and ALMAp^l ( down to sm> 
arcsecond resolution at submillimeter wavelengths 
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BLAST 


Source 


I 


b 
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$850 


S500 


$450 


$350 


$250 


$ioo 


$50 


ID 


Name 


(Degree) 


(Degree) 


(Jy) 

y^j ) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 

y^j ) 


(Jy) 


(Jy) 


CO 


J203057+394248 


78.6760 


0.1902 






22.0 ± 3.9 




85.8 ± 3.4 


110.7 ± 5.6 


15. 6 U 




CI 


J202925+401118 


78.8861 


0.7091 


8.1 


28.3 


151.5 ± 2.9 


226.9 


375.3 ± 8.0 


1029.3 ± 35.2 


5819.4 


5699.0 


C2 


J202904+401506 


78.8980 


0.8006 






90.7 ± 32.4 




79.1 ± 17.8 


59.9 ± 6.7 


40.3 


1.3 


C3 


J203329+394127 


78.9454 


-0.2135 






23.7 ± 1.3 




50.9 ± 2.9 


85.4 ± 5.0 


41. l u 


7.7 U 


C4 


J203323+394253 


78.9546 


-0.1860 






33.0 ± 1.3 




68.8 ± 2.9 


130.7 ± 5.6 


118. 2 U 


6.9 U 


C5 


J202653+403717 


78.9566 


1.3528 






5.1 ± 0.4 




7.7 ± 0.5 


11.3 ± 0.9 


11. 2 U 


6.1 u 


C6 


J203023+400926 


78.9702 


0.5401 






15.2 ± 0.9 




50.3 ± 2.4 


74.4 ± 4.4 




1.0 U 


C7 


J203112+400310 


78.9778 


0.3528 


6.6 


17.6 


133.8 ± 3.0 


187.6 


287.8 ± 5.1 


512.5 ± 10.1 


935.1 


688.1 


C8 


J202957+401533 


79.0029 


0.6686 


4.0 




42.0 ± 2.5 




80.6 ± 4.8 


150.8 ± 10.2 


876. 7 U 


253. U 


C9 


J203242+395409 


79.0262 


0.0331 






9.7 ± 1.3 




20.9 ± 2.1 


43.3 ± 4.2 


40. 2 U 


39. 8 U 


CIO 


J203002+401637 


79.0279 


0.6645 


4.0 




91.3 ± 5.6 




145.5 ± 7.5 


198.4 ± 13.0 


195. 5 U 


51. 8 U 


Cll 


J203230+395738 


79.0499 


0.0986 






16.3 ± 1.4 




37.0 ± 3.0 


70.2 ± 5.3 


89.3 


52.0 


C12 


J203220+400021 


79.0687 


0.1497 






11.9 ± 1.6 




26.3 ± 4.4 


71.9 ± 7.4 


201. 8 U 


42. U 


C13 


J202903+402817 


79.0746 


0.9314 






10.9 ± 2.3 




28.0 ± 5.7 


46.2 ± 10.1 


10. 8 U 


18. 8 U 


C14 


J202720+404234 


79.0776 


1.3350 






8.9 ± 1.0 




19.8 ± 2.6 


32.8 ± 2.9 


116.3 


53.2 


C15 


J202904+403148 


79.1237 


0.9637 






107.1 ± 6.9 




175.6 ± 7.6 


259.2 ± 9.4 


292.2 


122. 7 U 


C16 


J203442+394448 


79.1312 


-0.3696 






95.9 ± 5.0 




167.0 ± 4.0 


286.2 ± 6.0 


540.8 


85.4 


C17 


J203326+395612 


79.1370 


-0.0594 




1.3 


5.9 ± 1.1 


9.4 


10.6 ± 1.6 


29.2 ± 3.4 


112. l u 


20. 3 U 


C18 


J203042+401921 


79.1388 


0.5894 






11.2 ± 1.4 




22.3 ± 3.0 


32.4 ± 3.5 


18. 2 U 




C19 


J203108+401718 


79.1606 


0.5016 






13.0 ± 1.1 




29.7 ± 3.1 


76.0 ± 4.1 


452.1 


178.2 


C20 


J203229+400848 


79.1987 


0.2105 






19.3 ± 3.5 




38.7 ± 6.0 


91.9 ± 11.3 


428. 5 U 


137.3 U 


C21 


J202655+405517 


79.2044 


1.5223 






22.7 ± 2.5 




53.7 ± 3.4 


94.8 ± 8.5 


12. 6 U 


2.3 U 


C22 


J202923+403603 


79.2173 


0.9549 






32.7 ± 2.7 




82.1 ± 9.5 


92.9 ± 11.3 


27.8 U 


8.7 U 


C23 


J203116+402231 


79.2462 


0.5321 






31.3 ± 3.2 




99.8 ± 22.9 


92.7 ± 12.6 




0.1 u 


C24 


J203138+401939 


79.2477 


0.4492 






15.5 ± 3.4 




10.1 ± 6.4 


9.0 ± 4.6 




69. 8 U 


C25 


J203146+401843 


79.2499 


0.4198 




6.8 


67.4 ± 2.1 




133.5 ± 7.7 


121.1 ± 7.2 






C26 


J203158+401833 


79.2702 


0.3877 


3.2 


10.1 


44.2 ± 2.2 


135.1 


70.1 ± 4.1 


87.3 ± 4.5 


319. 5 U 


45. 7 U 


C27 


J202806+405127 


79.2822 


1.3037 






32.3 ± 1.8 




124.4 ± 8.3 


292.5 ± 18.7 


1312.2 


1217.3 


C28 


J203227+401554 


79.2903 


0.2855 


5.4 




44.3 ± 4.9 




92.2 ± 8.4 


307.5 ± 31.2 


2885.4 


2740.1 


C29 


J202738+405706 


79.3081 


1.4293 






8.0 ± 1.1 




12.8 ± 2.2 


26.2 ± 2.4 


129. 3 U 


21. 4 U 


C30 


J202809+405301 


79.3099 


1.3099 






91.9 ± 3.1 




199.9 ± 7.8 


453.2 ± 16.4 


2654.2 


2324.4 


C31 


J203518+395310 


79.3107 


-0.3767 






36.9 ± 2.1 




100.8 ± 5.2 


141.7 ± 7.1 


651. 7 U 


289. 6 U 


C32 


J203233+401651 


79.3145 


0.2796 


5.1 




107.2 ± 5.5 




275.0 ± 9.8 


666.5 ± 31.4 


8224.0 


6221.7 


C33 


J202818+405309 


79.3275 


1.2893 






21.7 ± 3.3 




22.4 ± 4.8 


84.9 ± 10.8 


3271. 4 U 


1948. l u 


C34 


J203525+395328 


79.3278 


-0.3913 






14.1 ± 1.9 




33.8 ± 3.2 


93.7 ± 7.1 


272. 3 U 


161. 9 U 


C35 


J202823+405234 


79.3292 


1.2704 






6.9 ± 11.2 




19.2 ± 12.7 


73.1 ± 34.6 


2234. 4 U 


913. 7 U 


C36 


J203222+401955 


79.3337 


0.3394 


6.6 


20.4 


123.8 ± 3.3 


172.8 


238.5 ± 5.5 


386.5 ± 13.3 


1139.2 


277.0 


C37 


J202802+405636 


79.3459 


1.3622 






11.6 ± 1.1 




19.0 ± 2.1 


42.8 ± 3.7 


382. 5 U 


65. 3 U 


C38 


J203349+400832 


79.3471 


0.0025 






29.8 ± 0.9 




68.3 ± 1.9 


172.1 ± 3.9 


978.1 


988.4 


C39 


J203530+395438 


79.3531 


-0.3928 






17.6 ± 1.7 




57.9 ± 5.0 


74.2 ± 8.7 


444. 6 U 


270. 6 U 


C40 


J203508+395725 


79.3492 


-0.3102 






32.0 ± 3.2 




44.3 ± 2.6 


103.2 ± 9.2 


23. 8 U 


1.3 U 


C41 


J202818+405638 


79.3746 


1.3231 






4.2 ± 0.4 




12.4 ± 2.6 


51.0 ± 5.1 


479. 9 U 


121. 4 U 


C42 


J203530+400053 


79.4371 


-0.3311 






21.4 ± 3.0 




32.8 ± 2.0 


44.8 ± 3.2 


101. 6 U 


11. 6 U 


C43 


J202952+404842 


79.4415 


1.0049 






21.2 ± 1.0 




50.9 ± 1.9 


81.5 ± 2.9 


245.8 


137.6 


C44 


J202754+410505 


79.4460 


1.4651 






13.5 ± 1.1 




39.4 ± 3.4 


82.3 ± 6.7 


100.4 


15.7 


C45 


J203616+395728 


79.4807 


-0.4836 






19.7 ± 3.3 




63.7 ± 3.5 


137.7 ± 5.7 


87.1 


7.1 


C46 


J203412+401634 


79.4985 


0.0231 






17.7 ± 1.4 




32.8 ± 2.0 


71.3 ± 8.5 


7.6 U 


2.7 U 


C47 


J203417+401840 


79.5360 


0.0314 






16.1 ± 0.8 




32.0 ± 1.6 


26.2 ± 2.2 


19. 4 U 


12. 4 U 


C48 


J203423+403039 


79.7073 


0.1358 






28.6 ± 3.3 




77.8 ± 10.7 


55.7 ± 4.7 


46. 4 U 


5.2 U 


C49 


J203311+404138 


79.7185 


0.4278 






6.4 ± 0.4 




11.8 ± 0.6 


29.2 ± 1.1 


102.4 


49.2 


C50 


J203610+401716 


79.7318 


-0.2685 






14.7 ± 1.1 




36.6 ± 2.4 


45.3 ± 3.4 






C51 


J203050+410224 


79.7340 


0.9909 






24.3 ± 0.8 




50.9 ± 1.7 


79.5 ± 2.3 


95.2 


55. 4 U 


C52 


J203119+410049 


79.7670 


0.9011 






6.0 ± 0.6 




13.4 ± 0.9 


22.2 ± 1.7 


33. 4 U 


19. 7 U 


C53 


J202958+411639 


79.8299 


1.2628 






10.1 ± 0.6 




34.6 ± 1.4 


53.6 ± 2.2 


133.0 


84.4 


C54 


J203134+410357 


79.8364 


0.8949 






8.4 ± 0.8 




16.3 ± 1.6 


14.8 ± 1.5 


53. 6 U 


1.6 U 


C55 


J203029+411556 


79.8764 


1.1790 






90.0 ± 2.8 




200.4 ± 3.5 


442.4 ± 7.0 


1470.4 


820.4 


C56 


J203208+410524 


79.9188 


0.8234 






90.0 ± 2.8 




200.4 ± 3.5 


442.4 ± 7.0 






C57 


J203148+410836 


79.9241 


0.9065 






-21.5 ± 26.8 




3.9 ± 1.4 


15.6 ± 3.0 


166. U 


79. 9 U 


C58 


J203319+405858 


79.9643 


0.5814 






4.3 ±0.9 




16.5 ±2.9 


25.6 ±4.1 


112.4 


51.6 


C59 


J203223+410748 


79.9787 


0.8099 






7.3 ±0.7 




24.5 ±0.9 


40.6 ± 1.4 


258.5 U 


242.7 U 


C60 


J203201+411053 


79.9789 


0.8963 






32.1 ±2.1 




72.6 ±5.3 


115.4 ±8.3 


579.5 U 


51. 7 U 


C61 


J203216+410860 


79.9813 


0.8397 






63.5 ±11.3 




133.0 ± 19.8 


206.6 ±35.5 


889.9 


184.1 


C62 


J203437+405147 


80.0149 


0.3127 






87.2 ±4.1 




171.6 ±7.4 


311.9 ± 13.1 


63.0 U 


70.2 U 


C63 


J203745+402529 


80.0234 


-0.4268 






3.9 ±0.6 




9.4 ± 1.5 


15.7 ± 1.6 


41. 7 U 




C64 


J203443+405308 


80.0457 


0.3090 






20.9 ±1.1 




61.7 ±5.0 


96.7 ± 10.9 


126.0 


59.7 


C65 


J203309+411603 


80.1752 


0.7749 






10.5 ±0.7 




17.1 ±0.9 


36.0 ± 1.5 


265.8 


70.0 


C66 


J203414+410808 


80.1903 


0.5331 






57.6 ±3.5 




94.2 ±3.4 


168.6 ±7.3 


70.8 


76.6 


C67 


J203345+411415 


80.2188 


0.6657 






7.9 ±1.1 




12.7 ±0.9 


22.8 ± 1.6 


208.1 


73.9 


C68 


J203520+410551 


80.2851 


0.3429 






16.0 ±0.7 




43.0 ±1.3 


79.9 ±2.1 


23.6 U 




C69 


J203112+414229 


80.3133 


1.3307 






11.4 ±0.7 




39.2 ±3.1 


45.2 ±3.4 


31.4 


4.0 


C70 


J203537+410615 


80.3220 


0.3050 






5.1 ±0.4 




13.0 ±1.1 


18.0 ± 1.1 


32.3 


60.7 U 
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TABLE 1 - Continued 
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Qn7 c 
oU / .5 


OO Q 

oz.o 


C167 


TonQon/i i /lOQO/ic 
J zUooU4+4zoy45 


oi oQnn 

oi.ooyy 


0.8813 


5.4 




A C T _L O 7 

45.7 zt z.7 




1 Qn _L ^ c 

loy.z ± b.5 


OI C O _L 1 1 Q 

zl5.z zt ll.o 


812. 8 U 


99. U 


C168 


Tor\97f;7 i /io/incn 
JzUo757+4z4U5U 


o i o a i n 
ol.o41U 


n onno 

u.yuyo 






Q Q 1 _L Q Q 
OO.l ± 0.0 




on /i n _L 1 n 
z(J4.U zt lU.z 


Q A 7 n _L 1/11 

o47.y ± 14.1 


404.6 


26.5 


r^i pci 


TOOQQQ7 1 /IOQ7QO 

J zUooo / +4zo / oo 


OI Q7QQ 


u. / / yo 


ob.o 


Q7 Q 
O / .O 


00 k J_ q n 
zz.5 zt o.U 


1 nno n 
lUUo.U 


QQ /? 1 /? 

00. b zt b.o 


/in c j_ 1 1 
4y.5 zt lz.1 


1 C /I Q Q O 

15400.Z 


1 1 n7n n 
iiy ( U.U 


Ul / u 


J ZUob44H-4z55U / 


OI 

ol.oybl 


l.zoUy 






c n J_ n 
515. y zt y.o 




1 1 /i n ^ j_ 00 
114U.b zt zo.z 


Qono A 1 A P 1 
ozUz.4 zt 4b. 1 


1 OK 

1Z5.5 


1 Q ^u 

lo.b 


f^l 71 


JzU4zo < H-4zzblo 


OO 1 7Q /I 

oz.!7o4 


n nonn 
U.UoUy 






10.6 zt 1.5 




Q O O _L Q 7 
OO.O ± 0.7 


on /i _i_ q n 

oU.4 ± o.y 


219.5 


154.5 


L/17z 


TOO/lOQ/l I /I0070£? 

J zU4zo4H-4zz f zb 


oo 1 ooc 

oz.loob 


0.1005 






24.8 zb 1.6 




55.7 zt 3.3 


1 n7 7 _L /l A 

1U7.7 ± 4.4 


98. 9 U 


Q7 1 U 

o7.1 


/^1 7Q 


JzU4UloH-4zo4oy 


OO OOI A 

OZ.Z014 


n 700Q 

U. ( zzo 






oq n J_ 1 A 
zo.U zt 1.4 




a n a J_ 
4U.4 zt z.o 


re n 1 Q A 

55. U zt o.4 


1 nu 

oi.y 


A 7U 
4. i 


C174 


T9D4D1 7+49o694 

JZVJ^fcU-L 1 |^ 4: Z ; J VJ Z/ ^± 


so SD7S 


0.7280 






17.5 zb 0.7 




OC C 1 1 C 

. _l_ J. . (J 


OI O _L C O 
Z/ J. . O _l_ . Z 


126. l u 


40. 6 U 


C175 


J204227+424151 


82.3591 


0.2663 






14.8 zb 0.7 




31.7 ± 1.6 


56.6 zt 2.1 


22. 3 U 


3.7 U 


C176 


J204159+424812 


82.3909 


0.3972 






65.6 zb 5.6 




91.2 zb 5.3 


147.6 zb 8.1 


185.1 u 


41.4 U 


C177 


J204211+425258 


82.4752 


0.4182 






13.9 zb 1.1 




31.4 =b 2.3 


80.8 zb 6.1 


526.3 


80.5 


C178 


J204209+425519 


82.5020 


0.4479 






68.2 =b 4.7 




97.4 =b 6.6 


222.5 zb 12.8 


618.1 u 


77.6 U 


C179 


J204341+424355 


82.5271 


0.1077 






17.3 zb 1.7 




10.0 zb 5.3 


78.2 =b 9.9 


54.9 U 


4.1 u 


C180 


J204234+425649 


82.5689 


0.4029 






80.7 zb 15.1 




88.9 zb 5.6 


141.9 zb 10.1 


573.4 


769.8 


C181 


J204346+424707 


82.5787 


0.1285 






14.3 zb 1.4 




53.0 zb 4.1 


125.3 zb 6.6 


9.4 U 




C182 


J204222+425831 


82.5695 


0.4482 






31.4 zb 3.3 




27.3 zb 2.2 


32.1 zb 3.9 


500.2 


179.1 


C183 


J204329+424954 


82.5826 


0.1985 






20.5 zb 4.7 




103.8 zb 6.3 


222.2 zb 14.0 


168.4 


95.5 



u Upper limit. 
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TABLE 2 

Clusters and Stellar Groups 



Cluster a 


I 


5 


Other 


ID 


Degree 


Degree 


Names 


LK08 


79.302 


0.286 


DB10, ECX6-25 


LK09 


79.301 


1.291 


DB11, ECX6-18 


DB12 s s 


79.879 


1.178 


ECX6-21 


ECX6-27 


80.354 


0.728 




DR18 s s 


80.357 


0.450 




LK11 


80.935 


-0.167 


DB13, DB14 S § 


LK12 


81.298 


1.096 


DB15 


LK13 


81.445 


0.483 


DB16 


LK14 


81.441 


1.122 




DB17 


81.566 


-0.719 




LK15 


81.610 


0.142 




DB18 S § 


81.661 


-0.017 




DR21 


81.680 


0.540 




DB19 s s 


81.711 


0.582 




W75N 


81.870 


0.779 




DB20 


81.896 


0.797 




DB22 


82.567 


0.404 





a LIC DB, and ECX6 a re the cluster I Ds by 
|Le Duigou & Knodlseder (2002), Dutra fe Bica] 
(2U0l]>, andfUomeron et ai. (20U2), respectively. 
^ Stellar group. 



TABLE 3 

Structures observed at 1.2 mm associated with BLAST clumps 



BLAST 




Core a 




Clump a 


^core 


V , a 
^emmp 


ID 




ID 




ID 


g cm -2 


g cm -2 


CI 




S26 




C1-S5 


0.58 


0.33 


C7 




S30, S31 




C1-S7 


0.21 


0.18 


C8 




S29 




C1-S6 


0.08 


0.08 


C26 




S34 




C1-S9 


0.08 


0.08 


C28 




S39 






0.12 




C32 




S40, S41 




Cl-Sll 


0.13 


0.11 


C36 




S36, S37, S38 




C1-S10 


0.16 


0.12 


C85 




N5, N6 




C1-N3 


0.15 


0.15 


C87 




N14 




C1-N6 


0.48 


0.14 


C88 




N10 




C1-N4 


0.43 


0.18 


C90 




N17 






0.15 




C92 




N58, N59 




C1-N17 


0.22 


0.18 


C104 




N62 




C1-N20 


0.08 


0.09 


C107 




N63 




C1-N21 


0.96 


0.28 


C115 




Nl, N2 




Cl-Nl 


0.56 


0.29 


C116 




N70 






0.07 




C119 




N12, N13 




C1-N5 


0.58 


0.21 


C126 




N18 




C1-N9 


0.07 


0.07 


C127 




N29 




Cl-Nll 


0.06 


0.08 


C128 




N16 






0.13 




C132 




N67, N69 




C1-N23 


0.05 


0.09 


C134 




N26 






0.06 




C135 




N61 






0.07 




C139 




N8 






0.12 




C141 




N65, N64 




C1-N22 


0.14 


0.10 


C145 




N68 






0.29 




C149 




N72 






0.05 




C155 




N42, N46, N47 




Cl-15 


2.30 


1.67 


C157 


N36, 


N38, N41, N44, 


N48 


C1-N14 


1.92 


1.24 


C160 


N37, 


N43, N51, N53, 


N54 


C1-N16 


0.45 


0.45 


C163 




N20 






0.18 




C166 




N9 






0.05 




C167 




N21, N22, N24 




C1-N10 


0.08 


0.10 


C169 




N30, N31, N32 




C1-N13 


2.49 


1.27 



a Core, Clump ID, and surface density are from |Motte et al7| ( |2Q07| ) 



BLAST Survey of Cygnus X 
TABLE 4 
RESULTS FROM SED FITS 



BLAST 


T 


M 


s 


L 


^bol 


ID 


(K) 


(100 M ) 


(g cm -2 ) 


(100 L ) 


(100 L ) 


CO 


15.5 ±0.7 


1.8 ±0.4 


0.033 


1.4 ±0.2 




CI 


40.1 ±4.0 


1.9 ±0.4 


0.121 


277.2 ± 100.7 


489.4 


C3 


16.6 ± 1.1 


1.4 ±0.3 


0.021 


1.6 ±0.4 




C4 


17.5 ±0.7 


1.7±0.3 


0.022 


2.6 ±0.4 




C5 


13.2 ± 1.3 


0.5 ±0.2 


0.022 


0.1 ±0.0 




C6 


19.6 ±0.6 


0.7±0.1 


0.009 


1.8 ±0.2 


2.0 


C7 


25.8 ±0.8 


2.6 ±0.3 


0.084 


33.4 ±3.3 


50.2 


C8 


17.3 ± 1.6 


2.2 ±0.6 


0.038 


3.1 ± 1.0 




C9 


20.4 ±2.0 


0.4 ±0.1 


0.006 


1.3 ±0.4 




CIO 


15.7 ± 1.3 


4.1 ±1.0 


0.051 


3.4 ± 1.0 




Cll 


25.7 ±0.8 


0.4 ±0.1 


0.008 


4.4 ±0.4 




C12 


20.9 ±0.6 


0.4 ±0.1 


0.007 


1.8 ±0.2 




C13 


15.7 ±0.8 


0.8 ±0.2 


0.014 


0.7 ±0.1 




C14 


29.9 ± 1.7 


0.1 ±0.0 


0.003 


3.8 ±0.6 


4.5 


C15 


20.8 ±0.9 


2.7 ±0.5 


0.024 


10.5 ± 1.1 


15.3 


C16 


23.7 ±0.6 


2.1 ±0.3 


0.096 


16.5 ± 1.4 




C17 


21.5 ± 1.8 


0.2 ±0.0 


0.005 


1.0 ±0.3 




C18 


15.1 ±2.0 


0.8 ±0.3 


0.015 


0.5 ±0.2 




C19 


32.5 ± 1.1 


0.2 ±0.0 


0.003 


9.8 ±1.2 


10.1 


C20 


23.6 ±3.8 


0.5 ±0.2 


0.006 


4.1 ± 1.8 




C21 


15.2 ±0.6 


1.6 ±0.3 


0.009 


1.1 ±0.1 




C22 


14.7 ± 1.2 


2.5 ±0.7 


0.021 


1.5 ±0.3 


2.3 


C23 


15.0 ± 1.6 


2.2 ±0.8 


0.015 


1.5 ±0.5 




C25 


15.1 ± 1.2 


2.4 ±0.5 


0.018 


1.6 ±0.5 




C26 


13.4 ±0.9 


3.7 ±0.6 


0.059 


1.3 ±0.3 




C27 a 


40.1 ±4.2 


0.4 ±0.1 


0.007 


62.3 ±24.9 


55.4 


C28 


47.0 ±9.7 


0.4 ±0.1 


0.009 


144.3 ±211.5 


113.5 


C29 


17.3 ±2.7 


0.4 ±0.2 


0.007 


0.5 ±0.2 




C30 a 


37.7 ±4.6 


1.1 ±0.5 


0.020 


110.6 ±37.2 


102.5 


C31 


18.0 ±2.3 


1.8 ±0.6 


0.033 


3.3 ±1.2 




C32 


46.2 ±2.8 


1.1 ±0.2 


0.024 


330.7 ±72.1 




C34 


25.4 ± 1.3 


0.4 ±0.1 


0.006 


4.3 ±0.7 




C36 


25.5 ± 1.1 


2.6 ±0.4 


0.039 


30.6 ±4.3 




C37 


18.8 ±2.5 


0.4 ±0.2 


0.011 


1.0 ±0.4 




C38 


37.6 ± 1.9 


0.4 ±0.1 


0.014 


38.3 ±6.0 


48.7 


C39 


19.4 ±2.8 


0.8 ±0.3 


0.010 


2.0 ±1.3 




C40 


15.5 ± 1.2 


1.7 ±0.5 


0.009 


1.3 ±0.3 




C41 


35.1 ± 1.9 


0.1 ±0.0 


0.002 


4.4 ±0.9 


5.4 


C42 


13.0 ± 1.6 


2.0 ±0.8 


0.027 


0.6 ±0.2 




C43 


28.9 ± 1.1 


0.4 ±0.1 


0.010 


8.4 ±0.9 




C44 


23.0 ±0.5 


0.4 ±0.1 


0.003 


3.0 ±0.2 




C45 


18.4 ±0.3 


1.1 ±0.1 


0.012 


2.2 ±0.2 




C46 


14.9 ±0.7 


1.2 ±0.3 


0.004 


0.8 ±0.1 




C47 


11.0 ±0.9 


2.6 ±0.7 


0.032 


0.3 ±0.1 


0.9 


C48 


11.8 ± 1.2 


3.8 ±1.3 


0.025 


0.7 ±0.2 




C49 


27.0 ±0.9 


0.1 ±0.0 


0.005 


2.1 ±0.2 


2.8 


C50 


15.3 ± 1.7 


1.0 ±0.3 


0.010 


0.8 ±0.2 




C51 


22.2 ± 1.0 


0.7 ±0.1 


0.026 


3.7 ±0.5 


5.0 


C52 


17.8 ±3.0 


0.3 ±0.1 


0.005 


0.5 ±0.1 


1.2 


C53 


27.8 ± 1.1 


0.2 ±0.0 


0.005 


4.2 ±0.5 


4.9 


C54 


11.5 ± 1.1 


1.2 ±0.4 


0.055 


0.2 ±0.0 




C55 


30.2 ± 1.3 


1.6 ±0.2 


0.117 


48.2 ±6.2 


52.9 


C57 


28.8 ± 1.4 


0.1 ±0.0 


0.004 


2.3 ±0.4 


3.6 


C58 


28.3 ± 1.0 


0.2 ±0.0 


0.007 


3.3 ±0.4 


4.1 


C59 a 


17.6 ± 1.9 


1.7 ±0.5 


0.053 


2.6 ±0.9 


6.8 


C60 a 


16.4 ± 1.3 


3.8 ±1.6 


0.021 


4.0 ±2.4 




C61 a 


25.7 ±0.9 


1.8 ±0.2 


0.029 


22.8 ±2.9 


24.8 


C63 


17.5 ±0.7 


1.1 ±0.3 


0.006 


1.7 ±0.2 




C64 


30.6 ± 1.7 


0.1 ±0.0 


0.006 


4.3 ±0.6 


6.0 


C65 


24.4 ±0.9 


1.1 ±0.2 


0.008 


10.0 ±0.9 




C66 


30.5 ± 1.7 


0.1 ±0.0 


0.003 


2.6 ±0.4 


8.0 


C67 


25.5 ±0.8 


0.4 ±0.0 


0.009 


4.9 ±0.5 


5.6 


C68 


17.2 ± 1.5 


0.6 ±0.1 


0.006 


0.8 ±0.3 




C69 


22.4 ±0.6 


0.1 ±0.0 


0.008 


0.9 ±0.1 




C70 


34.3 ±2.3 


0.0 ±0.0 


0.001 


2.6 ±0.4 




C71 a 


30.7 ± 1.3 


0.9 ±0.1 


0.016 


30.6 ±4.5 


31.3 


C73 


33.4 ± 1.5 


1.6 ±0.2 


0.030 


84.3 ±17.1 


87.6 


C74 


23.2 ±2.6 


0.3 ±0.3 


0.006 


2.3 ±0.7 


3.1 


C75 


34.8 ±2.6 


0.1 ±0.0 


0.002 


6.2 ±1.3 


6.3 


C76 


28.4 ± 1.0 


0.8 ±0.1 


0.016 


18.3 ±2.1 


21.8 



Continued on Next Page. . . 
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TABLE 4 - Continued 



BLAST 


T 


M 






-^bol 


1JJ 




( 'i nn A/r ^ 
(^1UU IVLq J 


(g cm -2 ) 


/ 1 nn t \ 
(1UU L j 


I 'i nn t 
^1UU L/0 


070 

U78 


oi n _L o n 
Zl.9 zt z.9 


0.2 zb 0.1 


0.003 


0.8 zb 0.3 




U79 


07 /I _L 1 /I 

Z7.4 zt 1.4 


0.2 zb 0.0 


0.002 


2.8 zb 0.4 




U8U 


on n _i_ i n 
oU.U ± l.U 


1.5 zb 0.2 


0.056 


43.4 zb 5.1 


48.9 


Lv81 


13.5 zt 1.2 


0.9 zb 0.3 


0.017 


0.3 zb 0.1 




08z 


11.5 zb 1.3 


1.5 zb 0.8 


0.037 


0.2 zb 0.1 




U83 


28.8 d= 1.1 


0.5 zb 0.1 


0.007 


10.8 zb 1.3 


11.6 


Lv84 


OQ A _i_ n o 

Z3.4 ± U.8 


0.5 zb 0.1 


0.006 


o 1-7 1 n a 

6.1 ± U.4 




085 


on n I i n 

z9.9 ± l.U 


1.2 zb 0.1 


0.053 


35.4 zb 5.0 


51.8 


U8b 


16.0 zb 1.9 


2.3 zb 0.8 


0.025 


2.1 zb 0.6 




U87 


On. o _L i n 
ZD. 8 zt l.U 


1.8 zb 0.2 


0.058 


oo n _i_ a n 
Z8.U ± 4.9 


33.6 


nno 

O88 


32.5 zb 1.1 


1.9 zb 0.2 


0.091 


Of 7 1 1 1 O 

85.7 ± 11.3 


106.0 




25.7 zb 1.2 


0.1 zb 0.0 


0.002 


1.8 zb 0.3 






1 £J O _l_ O Q 

lb. 8 zt 2.6 


3.3 zb 1.7 


0.019 


4.0 zb 1.3 




(^yi 


17.2 zb 1.2 


1.5 zb 2.3 


0.008 


2.0 zb 0.4 




nna 


39.3 zb 5.5 


1.9 zb 0.7 


0.024 


o a n n I n i o 

z4U.b ± 91.8 


323.2 


nno 


on n _i_ o i 
ZZ.U zt 2.1 


1.8 zb 0.5 


0.015 


9.7 zb 2.1 




one: 


O A A _l_ 1 7 

z4.4 zt 1.7 


0.3 zb 0.1 


0.001 


2.4 zb 0.6 






21.0 zt 1.2 


2.6 zb 0.5 


0.013 


inn i on 

lU.b ± z.U 




nnn 


19.1 zb 0.5 


2.1 zb 0.3 


0.009 


5.1 zb 0.4 




f(nn 

uyy 


2b. 6 zt U.9 


1.2 zb 0.2 


0.007 


1 Q 7 _i_ o n 

13.7 zt z.U 




C100 


o o r inn' 

23.5 ± U.b 


0.3 zb 0.0 


0.002 


1.9 zb 0.2 




OlUl 


on o i n o 

zb.3 ± U.8 


1.6 zb 0.2 


0.013 


o i n I o a 

zl.9 zt z.4 


23.0 


C102 


OO A _l_ n 7 

zz.4 ± U.7 


0.6 zb 0.1 


0.007 


3.4 zt 0.3 




C104 


OO o i n n 

z8.z ± u.y 


1.8 zb 0.2 


0.027 


O 7 7 1 O O 

37.7 ± 3.8 


44.7 


ninr 

01U5 


19.9 zt 1.7 


0.7 zb 0.2 


0.007 


2.1 zb 0.6 




C106 


15.6 zb 1.3 


2.0 zb 0.5 


0.018 


1.6 zb 0.4 




C107 


oo o _i_ n a 
22.2 zt U.b 


1.3 zb 0.2 


0.071 


7.3 zb 0.7 




C108 


15.4 ± U.b 


6.4 zb 1.1 


0.029 


4.8 zb 0.7 




C109 


i a a _i_ o o 

Ib.b zt z.8 


0.7 zb 0.2 


0.003 


0.8 zb 0.4 




U11U 


oo 7 _L n n 
z8.7 zt U.9 


0.7 zb 0.1 


0.006 


15.5 zb 1.5 




cm 


23.6 zt 1.4 


0.6 zb 0.1 


0.008 


5.0 zb 0.9 




Ollz 


19.5 zb 0.5 


1.4 zb 0.2 


0.015 


3.9 zb 0.4 




Ullo 


17.7 zb 0.5 


1.1 zb 0.2 


0.012 


1.7 zb 0.2 




C114 


ion i no 

iy.b ± U.8 


1.0 zb 0.2 


0.021 


2.7 zb 0.4 




C115 


o c n I no 

zo.b ± U.8 


2.2 zb 0.3 


0.115 


27.5 zb 3.1 


31.4 


Ullb 


17.1 zb 0.9 


1.7 zb 0.2 


0.022 


2.2 zb 0.4 




U117 


oo i _L n n 

3z.i ± u.y 


0.8 zb 0.1 


0.005 


34. b ± 6.b 




0118 


1 O O 1 O A 

18.3 ± z.4 


0.2 zb 0.5 


0.006 


0.4 zb 0.2 




uny 


21.8 zb 0.5 


3.0 zb 0.5 


0.115 


15.2 zb 1.5 




C120 


17.1 zb 0.7 


1.3 zb 0.3 


0.009 


1.8 zb 0.3 




ni Q1 

Ulzl 


1 Q Q _l_ O O 

lo.8 zt 2.2 


1.3 zb 0.6 


0.033 


0.5 zb 0.3 




C122 


12.9 zt 1.2 


5.4 zb 1.5 


0.079 


1.5 zb 0.4 




C123 


oi n _L n n 

zl.U zt u.y 


1.4 zb 0.2 


0.008 


5.9 zb 0.8 




C124 


OO A _l_ 1 O 

zz.4 zt l.z 


0.2 zb 0.0 


0.001 


1.4 zb 0.3 




C125 


oo i I n o 

z8.1 ± U.8 


0.2 zb 0.0 


0.002 


4.4 zb 0.6 


5.0 


C126 


21.9 zt 1.1 


2.7 zb 0.4 


0.019 


1 A n _l_ Q A 

14. U zt 3.4 


15.1 


C127 


00.9 zt o.7 


1.7 zb 0.3 


0.006 


nr i _|_ 71 i 

95.1 zt 71.1 


91.3 


C128 


27.1 zt 1.6 


0.3 zb 0.1 


0.012 


5.4 zt 0.9 


6.9 


C129 


on o I o o 

zU.3 ± 2.6 


2.7 zb 0.7 


0.026 


9.1 zb 3.8 


25.7 


C130 


OQ O _L 1 O 

26.2 zt l.z 


0.5 zb 0.1 


0.008 


3.7 zb 0.7 




C132 


27.2 zt 1.5 


2.7 zb 0.5 


0.015 


a a n _i_ n o 

4b. U ± 9.3 




C133 


23.1 zb 1.3 


0.6 zb 0.2 


0.004 


4.5 zb 0.9 




C134 


OO T _L n 

zz.7 zt U.b 


1.0 zb 0.1 


0.009 


6.4 zb 0.6 




C135 


25.5 zt 1.3 


0.4 zb 0.1 


0.011 


5.3 zb 0.8 


7.5 


C136 


14.6 zb 1.1 


5.1 zb 1.2 


0.040 


2.9 zb 0.6 




C139 


on t _L n <~r 
20.7 zt U.7 


0.4 zb 0.1 


0.011 


1.4 zb 0.2 




C140 


oo n _L n ^ 
ZZ.U zt U.b 


1.7 zb 0.2 


0.009 


9.1 zb 1.0 




c~\ 1/11 
0141 


21.5 zt 0.6 


1.7 zb 0.2 


0.041 


8.3 zb 0.7 


8.9 


C142 


21.6 zb 1.1 


0.5 zb 0.1 


0.007 


2.3 zb 0.4 




C143 


14.6 zb 1.5 


1.6 zb 0.5 


0.034 


0.9 zb 0.3 




C144 


i n a I n t 

iy.4 ± U.7 


0.5 zb 0.1 


0.006 


1.4 zb 0.2 




0145 


21.9 zb 1.2 


1.1 zb 0.1 


0.026 


5.6 zb 1.2 




U14b 


13.8 zb 1.9 


5.8 zb 2.7 


0.027 


2.4 zb 0.8 




U14 f 


1 Q Q _|_ O C 

lo.o zt Z.O 


i <-7 j_ i n 

i. ( ± i.y 


n noo 
U.Uzz 


n a J_ n q 
U.b zt U.o 




C148 


21.9 zt 0.6 


0.4 =t 0.0 


0.002 


1.8 zt 0.2 




C149 


19.5 zb 0.8 


1.1 zb 0.2 


0.022 


3.0 zb 0.6 




C150 


19.3 zb 0.9 


0.3 zb 0.1 


0.004 


0.8 zb 0.1 




C151 


20.2 zb 0.9 


0.8 zb 0.2 


0.014 


2.7 zb 0.4 




C152 


18.9 zb 0.8 


0.6 zb 0.1 


0.006 


1.5 zb 0.2 




C153 


19.7 zb 2.5 


0.7 zb 0.3 


0.019 


2.2 zb 0.7 




C154 


17.4 zb 0.8 


1.8 zb 0.3 


0.013 


2.7 zb 0.4 




C155 


36.0 zb 4.4 


8.8 zb 2.4 


0.432 


697.6 zb 250.4 


776.6 


C156 


26.1 zb 0.7 


0.4 zb 0.1 


0.004 


5.5 zb 0.6 





Continued on Next Page. 
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TABLE 4 - Continued 



BLAST 


T 


M 


5] 




T i i 

-^bol 


liJ 




(^1UU Mq J 


(g cm -2 ) 


I'y nn t \ 
(100 L j 


( ~\ nn t \ 
(100 L J 


C157 


zo.7 ± z.U 


18.8 ± z.4 


1.019 


z91.o ± 108. o 


319.3 


C158 


z4.o ± U.8 


0.2 zb 0.0 


0.001 


1.5 zb 0.2 




C159 


Z5.0 zb U.8 


0.8 zb 0.1 


0.007 


17.0 zb 1.7 




C160 


o it r\ i o o 

zo.9 ± z.3 


11.2 zb 2.2 


0.903 


1 A CT O I A <~7 A 

145.3 ± 47.4 


152.0 


C161 


zl.4 zb 2.9 


0.4 zb 0.1 


0.004 


1.7 zb 0.7 




C162 


12 A zb 1.9 


2.4 zb 1.2 


0.035 


0.5 zb 0.2 




C163 


2b A ± 2.1 


0.3 zb 0.1 


0.010 


4.9 zb 2.3 


5.8 


C164 


ZO.O ± O.Z 


0.3 zb 0.1 


0.007 


4.1 zb 1.8 




C165 


27.3 zb 1.5 


0.8 zb 0.1 


0.012 


14.0 zb 2.5 


17.3 


C166 


24.0 ± O.o 


1.2 zb 0.2 


0.013 


10.0 zb 1.0 




01o7 


OTA O I O 1 

zU.z zb 2.1 


3.2 zb 0.8 


0.045 


i n a i o o 

10. o ± z.8 




C168 


24.4 zb 1.3 


0.4 zb 0.2 


0.014 


3.7 zb 0.5 


4.9 


C169 


OD.Z ± O.O 


7.6 zb 1.2 


0.689 


DZ4.4 ± zz9.4 


664.9 


C170 


zz.9 ± 0.7 


0.4 zb 0.0 


0.002 


2.5 zb 0.3 




C171 


o <i i inn 

27.1 ± 0.9 


0.5 zb 0.1 


0.011 


8.3 zb 1.0 


11.5 


C172 


13.7 zb 1.6 


2.0 zb 0.6 


0.069 


0.8 zb 0.3 




C173 


10.2 zb 1.1 


3.5 zb 1.2 


0.007 


0.3 zb 0.1 




C174 


18.6 zb 2.0 


0.7 zb 0.2 


0.039 


L4 =b 0.7 




C175 


13.3 zb 1.1 


5.5 zb 1.4 


0.034 


1.8 zb 0.5 




C176 


25.2 zb 1.7 


0.4 zb 0.1 


0.004 


4.0 zb 0.9 




C177 


24.6 zb 0.8 


1.3 zb 0.2 


0.013 


12.4 zb 1.5 




C179 


14.1 zb 2.1 


4.3 zb 2.3 


0.035 


2.0 zb 0.7 




C180 


42.6 zb 2.4 


0.2 zb 0.0 


0.006 


34.4 zb 6.5 


49.4 


C181 


9.3 zb 1.3 


6.2 zb 3.5 


0.084 


0.3 zb 0.1 




C182 


29.4 zb 1.0 


0.6 zb 0.1 


0.007 


15.3 zb 1.8 




C183 


25.0 zb 2.2 


0.7 zb 0.7 


0.023 


7.1 zb 1.0 


10.2 



a M and L not corrected for different distance than Cyg OB2 (see § |4.2| 
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TABLE 5 

Compact HII regions and protostars 



T2T A CT 
D LAo 1 


T 


IVl 




Source 


MAMdU 


Figure a 


ID 


(K) 


(100 M ) 


(100 L ) 


Name 


1.2 mm 




No 


CI 


40 


1.9 


489.0 


S26, AFGL2591 


UCH II 


7 




C7 


26 


2.6 


50.2 


S30, IRAS 20293+3952 


MIRQP 






C14 


30 


3.8 


4.5 


IRAS 20255+4032 




b 




C30 b 


38 


17.6 


1640.0 


IRAS 20264+4042 








C32 


46 


1.1 


330.7 


S41, IRAS 20306+4005 


HLIRPC 


7 




C36 


26 


2.6 


30.6 


S37, IRAS 20305+4010 


MIRQP 


7 




C38 


38 


0.4 


48.7 


IRAS 20319+3958 








C49 


27 


0.1 


2.8 










C55 


30 


1.6 


52.9 


IRAS 20286+4105 








C58 


28 


0.2 


4.1 


IRAS 20315+4046 








C64 


31 


0.1 


6.0 


IRAS 20328+4042 








C71 b 


31 


22.5 


782.5 


IRAS 20320+4115 








C73 


33 


1.6 


87.6 


IRAS 20333+4102 








C76 


28 


0.8 


21.8 


IRAS 20327+4120 








C80 


30 


1.5 


48.9 


IRAS 20332+4124 








C85 


30 


1.2 


51.8 


N6, IRAS 20343+4129 


HLIRPC 


10 




C87 


27 


1.8 


33.6 


N14, IRAS 20352+4124 


HLIRPC 


1U 




C88 


33 


1.9 


106.3 


N10, IRAS 20350+4126 


UCH II 


1U 




C92 


39 


1.9 


323.2 


N58, IRAS 20375+4109 


CH II 


11' 




C104 


28 


1.8 


44.7 


N62 


IRQP 


T2 




C107 


22 


1.3 


7.3 


N63 


MIRP 


T2 




C115 


26 


2.2 


31.4 


N3 


MIRQP 


y 






C119 


22 


3.0 


15.2 


N12 


MIRQP 


y 




C141 


22 


1.7 


8.9 


N65 


MIRQP 


i 


2 




C145 


22 


1.1 


10.7 


N68 


MIRQP 


rz 




C155 


36 


8.8 


776.6 


N46, DR21 


UCH II 


ii 




C157 


27 


18.8 


319.3 


N44, DR21(OH) 


MIRQP 


ii 




C167 


20 


3.2 


10.6 


N24 


IRQP 


ii 




C169 


36 


7.6 


664.1 


N30, W75N 


HCH II 


ii 





a Figure number of the images containi ng in dividual sources. 
b M and L corrected for distance (see § 4.2). 
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APPENDIX 
MULTI- WAVELENGTH PHOTOMETRY 



Flux densities at many supplementary wavelengths are presented in Table [T] and Table [6| They were obtained as 
follows. 

MAMBO 1200 jam. The MAMBO camera mounted on the IRAM 30-m telescope imaged all regions having column 



density greater than (A v > 15) at wavelength of 1.2 mm ( Motte et aL| 2QQ7| >. With its 11" beam, MAMBO has 
the ability to resolve multiple cores (if present) within BLAST clumps. Combining our data with MAMBO helps in 
mult i- wavelength photometry, especially when any cold source is undetected in the mid to far-infrared spectrum. To 
obtain consistent photometry, we first convolved the MAMBO map to 1', about the resolution of the BLAST maps. 
This has two effects: blending together any MAMBO cores; and bringing in flux from the more extended emission 
(halo or plateau) usually associated with these multi-scale structures. Gaussian photometry was performed on the 
convolved map. 

SCUBA 850 and 450 /am. The archival dat a from the Submillimeter Common User Bolometer (SCUBA) on the 
15-m James Clerk Maxwell Telescope (JCMT),|Di Francesco et al. (2008) sample many high extinction regions of the 

^ 1 , • T-H 11- ^ M ^ T .,1 1 . , J I I •- 1—1 -J J B-fl 1 — — II ^ - ^ 1 A - ^ , . 1 



Galactic Plane, including Cyg X, with diffraction-limited beams of 14.5 and 7.5", at 850 and 450 //m, respectively. 
Due to the higher spatial resolution, SCUBA can also identify numerous small-scale cores. The archival images consist 
of discontinuous smaller maps and so could not be co nvolved to BLAST resolut ion prior to photometry. Instead, we 



used flux densities using the SCUBA legacy catalogs of Di Francesco et al. (2008). We obtained flux densities from the 
Extended Map Object Catalog (EMOC) and added up all liux densities for sources within 1' of the BLAST source. 
Because of the small SCUBA maps, some BLAST sources are not completely scanned near the map edges and the 
resulting flux density is low. 

IRAS 100 and 60 jam. Bright dust emission is observed at 100 and 60 jam throughout Cyg X, and for most of the 
more luminous, hotter BLAST sources these IRAS bands cover the peak of the SED. The IRAS point source catalog 
does not have any entries for about half of our Cyg X region. This arose because of the very bright source DR21 
(C. Beichman, private communication). Nevertheless, the IRAS images that have been made subsequently are not 
seriou sly affected. We extr act flux densities at 100 and 60 jam from the Infrared Galaxy Atlas images (IGA; |Cao et aL 



|1997[ ), made using HIRES ( |Aumann et al.|1990| ), a resolution enhancing algorithm. There is some elongation ol sources 
across the dominant scan direction. The resolution at 100 and 60 jam is about 2' and 1', respectively, substantially 
better resolution than in the original ISSA product. Most of the isolated BLAST sources are warm enough to have 
a counterpart in the IGA map. BLAST counterparts at 100 jam were measured with a circular aperture of 2A' . 
Similarly, at 60 jam elliptical apertures of dimension 1.8' x 1.2' were placed with the major axis aligned across the 
scan direction. In the crowded regions, multiple Gaussian photometry (elliptical in the case of 60 jam) was carried out 
by fixing centroids of nearby sources. 

The counterparts to BLAST sources in the 60 jam images are often affected by extended cirrus-like structure and 
sometimes by deconvolution artefacts. Such sources show inconsistent structure between the 100 and 60 jam images. 
We preferred not to use such photometric data as detections to c onstrain the Wien part of the SED, rather using them 
as upper limits through a penalty function ( Chapin et al.|[2008| . In crowded regions, the flux density measurements 



of fainter sources are highly uncertain and for such cases we also used the photometric data as an upper limit. These 
upper limits are indicated by a superscript V in Table [T] 

MIPS 70 jam. Spitzer MIPS observations at 70 /am have better resolution and sensitivity than IRAS IGA 60 jam 
data and so are particularly useful for detecting or placing upper limits on colder and fainter BLAST sources. The 
MIPS im ages were destrip ed, which is advantageous for probing faint sources, and corrected at the pixel level for non- 
linearity ( Dale et al.|2007 ). Measurements using a Gaussian model of fixed FWHM 40" (typical extent of actual MIPS 



sources) at the positions of 39 faint BLAST sources where we have found no apparent counterparts in the MIPS image 
give an rms flux density of 5.7 Jy. This provides a very useful upper limit for these sources, constraining the SED on 
the Wien side of the peak. For flux densities of somewhat brighter sources, we first convolved the 70 jam map to 1', to 
be comparable to BLAST, and then made measurements using multi- Gaussian photometry with a linear background 
model. This provided useful detections for 56 sources and upper limits for the others. For six BLAST sources, pixels 
at the peaks of the bright MIPS counterparts are blanked; however, the IRAS data for these are sufficient to constrain 
the SED. The resolution of our convolved maps is also closer to the IGA 60 jam resolution, facilitating a comparison 
with isolated sources with a good IRAS detection. We find that the MIPS flux densities are slightly higher than at 
60 /am (Fqo/Fjo = 0.78), as expected for these SEDs. At the native resolutio n, without beam matching, the MIPS 
flux densities are systematically lower, as is discussed by Mottram et al.| ( 2010[) . 

MSX 8, 12, lj, and 21 jam. MSX ( [Price et al.|l2001[ ) observed in A^UTD, and E bands at 8.3, 12.1, 14.7 and 
21.3 /im, respectively. Among the four bands, A is the most sensitive and is dominated by diffuse emission from PAH 
molecules. It has angular resolution 18.3". 

MSX sources have relatively small positional uncertainties (within 4" - 5"). MSX sources were identified within a 
search radius of 30" about the BLAST coordinate. After visual inspection, we rejected some MSX "matches" that 
were either in a PDR and/or excited by some nearby OB stars, and so had nothing to do with the star formation 
history of the subject BLAST source. 

We extracted flux densities from the map in its native resolution directly by fitting multiple Gaussians, assuming 
a constant (not tilted) background. We find that our flux densities in all the bands are systematically higher than 
reported in the MSX point source catalog, and the deviation is more noticeable for less bright sources. This trend was 
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TABLE 6 

MSX COUNTERPARTS 



BLAST 


Aa a 


AS a 


s 8 


#12 


#14 


#25 


ID 


\ ) 


\ ) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


CI 


5.9 


3.0 


313.8 


568.0 


798.5 


1023.4 


C6 


-1.5 


6.0 


0.2 


0.8 


0.6 


1.6 


C7 


-7.8 


13.7 


17.1 


23.5 


23.4 


100.0 


C14 


11.6 


-5.1 


0.4 


1.0 


1.2 


2.4 


C15 


17.2 


-6.2 


4.6 


5.8 


3.7 


13.6 


C19 


4.4 


2.3 


1.8 


1.7 


0.9 


1.3 


C22 


-27.6 


15.4 


0.7 


1.3 


2.0 


2.5 


C27 


15.4 


-6.3 


2.1 


4.7 


3.0 


8.7 


C28 


34.4 


-12.6 


4.8 


9.1 


6.6 


37.7 


C30 


-35.0 


-12.5 


0.5 


4.0 


7.6 


13.8 


C38 


8.7 


3.6 


12.5 


13.4 


11.1 


64.3 


C41 


8.0 


-5.8 


1.0 


1.0 


0.7 


2.0 


C47 


-4.4 


-10.7 


0.2 


0.6 


1.3 


1.5 


C49 


3.9 


-0.7 


0.7 


1.4 


1.3 


2.7 


C51 


6.9 


-1.5 


0.2 


0.6 


0.7 


2.3 


C52 


11.6 


4.7 


0.4 


1.0 


0.4 


1.5 


C53 


4.1 


-4.8 


0.8 


1.6 


1.0 


1.4 


C55 


-18.3 


-6.0 


6.3 


8.1 


4.0 


13.2 


C57 


10.2 


6.4 


0.5 


0.8 


0.4 


1.8 


C58 


-3.7 


-6.3 


0.7 


1.2 


1.1 


2.5 


C59 


-5.3 


18.7 


0.7 


1.2 


2.3 


2.7 


C61 


17.2 


0.7 


2.7 


3.5 


1.8 


3.7 


C64 


2.1 


0.6 


2.6 


4.7 


4.2 


5.7 


C66 


2.1 


5.9 


9.1 


15.9 


20.2 


26.8 


C67 


-2.0 


-1.4 


1.0 


1.6 


0.4 


1.3 


C71 


0.3 


10.5 


1.0 


1.9 


1.4 


2.8 


C73 


-17.3 


-2.2 


3.2 


5.8 


7.0 


19.5 


C74 


28.1 


-12.9 


0.4 


0.8 


0.5 


1.3 


C75 


1.0 


3.1 


0.5 


0.6 


0.5 


1.2 


C76 


0.1 


-7.6 


7.9 


9.1 


8.2 


7.4 


C80 


6.9 


-5.2 


4.7 


5.7 


3.9 


18.5 


C83 


-1.2 


-6.4 


1.3 


1.6 


0.7 


1.5 


C85 


2.7 


0.8 


10.2 


17.1 


23.7 


90.9 


C87 


0.8 


9.6 


3.4 


6.3 


9.6 


19.7 


C88 


9.2 


-0.0 


5.4 


14.6 


26.2 


99.3 


C92 


-21.6 


0.2 


27.1 


66.1 


98.5 


634.1 


C101 


1.6 


-13.0 


1.5 


1.7 


1.6 


1.4 


C104 


-6.0 


12.4 


4.9 


5.7 


4.4 


26.6 


C115 


4.9 


8.8 


1.5 


2.7 


3.0 


3.7 


C125 


-1.2 


15.0 


0.6 


0.9 


0.7 


1.2 


C126 


-1.0 


32.9 


1.4 


1.7 


1.5 


2.0 


C127 


-13.5 


22.0 


1.9 


3.4 


2.3 


5.8 


C128 


-6.5 


16.7 


3.1 


4.0 


1.8 


3.9 


C129 


11.6 


18.3 


0.8 


1.0 


0.7 


1.3 


C135 


8.0 


-7.8 


0.8 


2.7 


5.2 


11.6 


C141 


2.1 


-13.7 


0.3 


0.6 


0.7 


1.4 


C155 


3.3 


15.6 


8.8 


29.6 


68.4 


543.9 


C157 


-31.6 


26.6 


9.2 


16.3 


21.6 


72.7 


C160 


62.8 


13.0 


2.8 


3.3 


3.4 


10.8 


C163 


-14.0 


16.5 


0.3 


0.8 


1.2 


2.1 


C165 


2.6 


31.5 


6.2 


6.1 


4.3 


10.4 


C168 


-10.4 


11.2 


0.5 


1.3 


2.4 


3.7 


C169 


19.9 


8.7 


15.0 


27.1 


55.4 


331.8 


C171 


-0.4 


-6.4 


3.2 


4.1 


3.3 


17.0 


C180 


0.9 


2.5 


11.5 


17.0 


23.8 


124.5 


C183 


3.2 


11.8 


2.9 


6.5 


10.0 


13.3 



a Offsets of MSX sources with respect to BLAST counter- 
parts. 



also observed by iMolinari et al. (2008) 



Results of MSX photometry are given in Table [6] In our present analysis, MSX flux densities in no way constrain 
temperature and total mass. We use MSX flux densities in the SEDs only to calculate a small correction to the 
bolometric luminosity by connecting data points piecewise continuously ( |Chapin et al.|[2008| . See for example, the 
SED for W75N in Figure [14{ note the "rising spectrum" from 4 to 21 /im, characteristic ol a MYSO. We find that this 
correction does not change the results significantly as the Cyg X sources detected by BLAST are quite cool. 
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